AD-A205  540 


iiu  copy 


CONTRACT  REPORT  BRL-CR-605 


BRL 


1938  -  Serving  the  Army  for  Fifty  Year*  -  1988 


DYNAMICS  OF  A  BALLOTING  PROJECTILE 
IN  A  MOVING  GUN  TUBE 


K.  A.  ANSARI 
J.  W.  BAUGH,  JR. 


DECEMBER  1988 


DTIC 

E!  C,  I  b.  Ef* 

Nfo  mar  2  3  1988  0 

%  I 

H 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 


U.S.  ARMY  LABORATORY  COMMAND 

BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MARYLAND 

8  9  ° 


T  C>  9 

1.  O  /V 


DESTRUCTION  NOTICE 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the 
originator. 


Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical 
Information  Service,  U.S.  Department  of  Carmerce,  Springfield,  VA  22161. 


The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department 
of  the  Army  position,  unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers'  names  in  this  report  does  not  con¬ 
stitute  indorsement  of  any  ccmnercial  product. 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-01 BB 


ir*rni.LiJH«ir. 

gjg 

Lfttl 


ITY  CLASSlFtCATIO 
UNCLASSIFIED 


RITY  CLASSIFICATION  AUTHORITY 


SlFlCATION/ DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBERS) 

QRL-CR-bOS 


1b  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

APPROVED  FOR  PUBLIC  RELEASE *, 
DISTRIBUTION  UNLIMITED • 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Battelle  Pacific  NW  Lab 


6c.  AOORESS  (Crty,  Sfafa,  »nd  ZIP  Cod*) 


6b.  OFFICE  SYMBOL  I  7a.  NAME  OF  MONITORING  ORGANIZATION 
(if  applicable)  I 


SLCBR-IB-M 


P.0.  Box  999 
Richland,  WA  99352 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 
MSB,  IBD,  BRL 


Be.  ADDRESS  (Crty,  Stete.end  ZIP  Cod*) 

Director 

U.S.  Army  Ballistic  Research  Lab 
Aberdeen  Proving  Ground,  MD  21005-5066 


1 1 .  title  (Includ*  Security  Classification) 


7b.  ADDRESS  (C/ty.  State,  and  ZIP  Cod*) 


8b  OFFICE  SYMBOL  I  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(if  applicable)  f 


10  SOURCE  OF  FUNDING  NUMBERS 

PROJECT  fTASK_ 

18aAB626  N° 


Dynamics  of  a  Balloting  Projectile  in  a  Moving  Gun  Tube 


12.  PERSONAL  AUTHOR(S) 

K.A.  Ansari  and  J.W.  Baugh,  Jr. 


13a.  TYPE  OF  REPORT 

CR 


13b.  TIME  COVERED 
FROM _ TO 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  IIS  PAGE  COUNT 
January  1986  I 


COSATI  COOES 


GROUP 


18  SUBJECT  TERMS  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 


I9y ABSTRACT  ( Continue  on  reverse  if  necessary  end  Identify  by  block  number) 

In  recent  years,  there  have  been  significant  strides  made  in  modeling  gun  dynamics. 

Current  capabilities  include  the  use  of  flexible  gun  tube  and  projectile  descriptions. 
However,  modeling  the  interface  conditions  between  the  tube  and  projectile  is  still  an 
area  of  very  active  research.  Highly  detailed  finite  element  models,  which  are  currently 
under  development  at  the  BRL  and  at  other  agencies  will  lead  to  a  much  better  understanding 
of  the  important  parameters  involved  in  the  balloting  of  projectiles.  These  models,  however, 
require  both  substantial  computer  resources  and  time.  Thus  there  is  a  need  for  the 
development  of  simplified  tube/projectile  interface  descriptions  which  both  realistically 
represent  the  interface  and  are  easily  implemented.  The  bourrelet/gun-bore  and  obturator 
models  incorporated  in  this  study  reasonably  satisfy  these  requirements. 

continued  on  next  page  .. 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
UUNU'SSIFIEO/UNLIMITSP  □  SAME  AS  RPT 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Edward  M.  Patton 


DO  Form  1473,  JUN  86 


21  ABSTRACT  SECURITY  CLASSIFICAT.ON 
□  OTIC  USERS  UNCLASSIFIED 


22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

(301)  278-6805  SLCBR-IB-M 


Previous  editions  are  obsolete 

% 


SECURITY  CLASSIFICATION  Of  THIS  PAC 


19.  ABSTUCT  -  coat'd 


a  nonlinear  dynamic  analysis  of  a  balloting'  projectile  is  presented  in  this  report.  The 
mathematical  model  used  is  a  six-degree-of-freedom  coupled  model  of  the  projectile  and  gun 
tube  system.  The  effects  of  obturator  flexibility  and  projectile  impact  with  the  gun  bore 
at  the  bourrelet  are  included  in  the  analysis,  and  the  nonlinear  differential  equations  of 
motion  of  the  system  are  derived  using  a  Lagrangian  formulation.  Gaussian  elimination  and 
'Nevmark's  cons tant-average-aecelerat ion  method  are  employed  to  obtain  a  solution.  Finally, 
numerical  results  for  some  specific  test  cases  are  obtained  and  discussed. 


ULs  ■ 


DYNAMICS  OF  A  BALLOTING  PROJECTILE  IN 
A  MOVING  GUN  TUBE 


by 

KA  Ansarl 
JW  Baugh,  Jr. 


Pacific  Northwest  Lab 
Richland,  Washington 


Prepared  for 

U.S.  Army  Ballistic  Research  Laboratory 
Aberdeen  Proving  Grounds 
Aberdeen,  Maryland 


January  1986 


TABLE  OF  GOT  BITS 

PAGE 

LIST  OF  ILLUSTRATIONS . S 

I.  INTRODUCTION . 7 

II.  LITERATURE  REVIEW . 9 

III.  SYSTEM  DESCRIPTION . 10 

IV.  TUB  SYSTEM  KINETIC  ENERGY . ....10 

V.  BOU REELS!/ GUN- BORE  CONTACT  MODEL . 13 

VI.  OBTURATOR  MODEL . 14 

VII.  TUB  SYSTEM  POTENTIAL  ENERGY . IS 

VIII.  STATIONS  OF  MOTION . IS 

IX.  SOLUTION  TECHNIQUE . 19 

X.  COMPUTER  CODE  DESCRIPTION . 21 

XI.  NUMERICAL  RESULTS  AND  DISCUSSION . 21 

XII.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY . 32 

ACKNOWLEDGEMENT . 34 

REFERENCES . 35 

APPENDIX  A . 37 

APPENDIX  S . 41 

APPENDIX  C . 45 

APPENDIX  D . Si 

LIST  OP  SYMBOLS . 77 

DISTRIBUTION  LIST . 79 


3 


LIST  OF  ILLUSTRATIONS 

Figure 

1  Displaced  Position  of  Dun  Bore . 11 

2  Projectile  with  Bourrelet  and  Obturator  within  Gun  Bore  .  11 

3  Body-fixed  Axes  and  Projectile  Displacements . 12 

4  Obturator  and  Bourrelet/Gun-Bore  Contact  Models  .  13 

5  Time  History  of  Gas  Pressure  on  Projectile . 23 

6  Projectile  Foundation  Moment . 23 

7  Time  History  of  Horizontal  Gun  Displacement  (X^) . 24 

8  Time  History  of  Vertical  Gun  Displacement  (Yg) . 24 

9  Time  History  of  X-Dlsplacement  of  Projectile  (xp)  .  25 

10  Time  History  of  Gun  Angular  Displacement  (e)  for  Case  la . 26 

11  Time  History  of  Gun  Angular  Displacement  (0)  for  Case  lb .  26 

12  Time  History  of  Gun  Angular  Displacement  fe)  for  Case  2a . 27 

13  Time  History  of  Gun  Angular  Displacement  (e)  for  Case  2b . 27 


14  Time  History  of  Y-Dlsplacement  of  Projectile  (yp)  for  Case  la  ...  28 

15  Time  History  of  Y -Displacement  of  Projectile  (yp)  for  Case  lb  ...  28 

16  Time  History  of  Y-Dlsplacement  of  Projectile  (yp)  for  Case  2a  ...  29 

17  Time  History  of  Y-Dlsplacement  of  Projectile  (yp)  for  Case  2b  ...  29 

18  Time  History  of  Yawing  Motion  of  Projectile  (a)  for  Case  la.  ...  30 

19  Time  History  of  Yawing  Motion  of  Projectile  (o)  for  Case  lb  ....  30 

20  Time  History  of  Yawing  Motion  of  Projectile  (a)  for  Case  2a  ....  31 

21  Time  History  of  Yawing  Motion  of  Projectile  (<*)  fur  Case  2b  ....  31 

22  Bourrelet  Impacting  Bore . . 44 


'j 


I.  INTRODUCTION 


The  yawing  or  wobbling  motion  of  a  projectile  within  a  gun  tube  Is  an 
Important  consideration  In  Internal  ballistics  and  Is  known  as  balloting. 

This  motion  Is  a  function  of  a  number  of  small,  dlfflcult-to-measure  parame¬ 
ters  such  as  manufacturing  tolerances,  lack  of  concentricity  of  the  engraving 
of  the  obturator,  projectile  and  tube  deformation,  obturation  of  the  propel¬ 
lant  gases  and  obturator  wear,  Certain  combinations  of  these  parameters  can 
potentially  cause  malfunctions  to  occur  which  include  poor  projectile  launch, 
damage  to  the  fuse  mechanism  and  explosion  In  the  bore.  The  yawing  motion  of 
the  projectile  leads  to  yaw  at  shot  exit,  and  the  resultant  Inaccuracy  of  the 
round.  It  can  also  lead  to  muzzle  wear  which  can  be  significant  for  large 
bore  diameters  and  muzzle  velocities.  As  a  wobbling  projectile  moves  down  a 
gun  barrel,  several  forces  act  on  It  Influencing  Its  motion.  These  can  be  due 
to  imbalances,  asymmetries,  clearances  and  deformation  In  the  projectile 
Itself,  engraving  and  frictional  forces,  propellant  gas  blowby,  tube  vibra¬ 
tion,  recoil,  and  coupling  between  the  projectile  and  the  tube  motion. 
Depending  upon  the  magnitudes  Involved,  these  can  lead  to  excessive  local  tube 
wear,  excessive  projectile  engraving,  malfunctioning  of  projectfle  components^ 
and  undesirable  in-bore  motion.  They  can  also  cause  Improper  exterior  ballis¬ 
tic  free  flight  of  the  projectile.  Differences  occur  between  the  point  of  aim 
and  the  point  of  Impact,  which  can  be  ascribed  to  the  motion  of  the  gun 
Itself,  the  aerodynamic  forces  acting  and  the  yawing  motion  of  the  projectile 
In  the  gun  bore.  Although  these  differences  may  be  only  as  large  as  a  few 
mils.,  they  must,  obviously,  be  minimized.  In  order  for  this  to  be  effected,  a 
realistic  prediction  of  the  yawing  motion  of  the  projectile  1r  the  gun  tube 
becomes  necessary. 
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II.  LITERATURE  SURVEY 


In  orderj to  predict  the  first  maximum  yaw  exterior  to  the  gun,  Reno1 
carries  out  a;  rigorous  Lagranglan  treatment  of  the  angular  motion  of  the 
projectile  In  the  bore.  He  assumes  that  the  projectile  Is  centered  at  the 
rotating  band,  the  plane  of  yaw  rotates  with  the  rifling,  the  Impact  of  the 
bourrelet  Is  normal  to  the  tube  and  rebound  Is  described  by  a  coefficient  of 
restitution.  Friction  between  the  projectile  and  the  bore  Is  Ignored.  The 
projectile  Is  approximated  by  a  single  degree  of  freedom  pendulum  system,  and 
a  closed-form  solution  Is  obtained.  However,  the  Initial  yaw  observed  during 
development  of  a  36-Inch  mortar  In  1944  was  greater  than  that  predicted  by 
Reno's  theory.  Furthermore,  scratches  on  the  bore  surface  Indicated  that  the 

projectile  precessed  In  a  direction  opposite  to  the  spin  Imparted  by  the 

2 

rifling.  In  an  attempt  to  provide  an  explanation,  Thomas  generalizes  Reno's 
approach  by  removing  the  constraint  on  the  orientation  of  yaw  and  deduces  the 
motion  of  the  plane  of  yaw.  His  results  show  little  difference  In  the  yaw; 
the  orientation  of  the  plane  of  yaw,  however.  Is  quite  different.  Darpas 
studies  a  situation  In  which  the  bourrelet  touches  the  bore  throughout  the 
projectile  travel  within  the  gun  tube  and  concludes  that  the  gyroscopic  effect 
due  to  the  projectile  being  cocked  In  the  bore,  Is  a  predominant  one.  The 
eccentricity  of  the  center  of  gravity  and  the  droop  of  the  tube  may  also  not 
be  negligible.  Perdreavll 1e4  presents  an  analysis  of  the  projectile  balloting 
problem  using  Lagrange's  equations  and  Euler  angles.  The  projectile  Is 
assumed  to  have  no  center  of  mass  offset  and  to  be  In  dynamic  balance  about 
the  principal  geometric  axes.  The  equations  are  written  to  represent  complete 


.4 


^F.V.  Reno,  "The  Motion  of  the  Axis  of  a  Spinning  Shell  Inside  the  Bore  of  a 
Gun,"  BRL  Report  Bo.  BRL-R-320,  1943  (AD#  491839). 

n 

L.E.  Thomas ,  "The  Motion  of  the  Axis  of  a  Spinning  Shell  Inside  the  Bore  of  a 
Gun,"  BRL  Report  No.  BRL-R -544,  1945  ( PB #  22102). 

^J.G.  Darpas ,  Translated  by  B.P.  Hitohcook,  "Transverse  Forces  on  Projectiles 
Which  Rotate  in  the  Barrel,"  BRL  Report  No.  BRL-MR-1208 ,  1959  (Memorial  de 
l ' artilleri.e  francaise ,  31:19,  No.  1,  1957)  (AD#  218873). 

4 

F.J.  Perareaville,  "Analysis  of  the  Lateral  Motion  of  a  Projectile  in  the  Gun 
Tube,"  Sandia  Laboratories,  Albuquerqi/e,  NM,  SC-RR-710071,  1971. 
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lateral  dynamic  freedom  of  the  projectile  bourrelet  within  the  bore.  Impact 
Is  Includeu  and  the  equations  are  set-up  so  that  the  motion  of  the  gun  tube 
can  also  be  Introduced.  Chu  and  Soechtlng^  present  the  same  theory  although 
they  resort  to  Euler's  dynamical  equations.  Perdreavllle^  extends  his  analy¬ 
sis  of  Reference  4  to  Include  the  effects  due  to  an  unbalanced  projectile  In  a 
rigid  gun  tube.  In  a  subsequent  document7  he  develops  equations  of  motion 
describing  the  lateral  motion  of  an  artillery  projectile  as  It  moves  down  an 
elastic  gun  tube  that  vibrates  laterally.  However,  he  does  not  generate  any 

O 

numerical  results  in  either  of  these  reports.  Langhaar  and  Boresl  analyze 
the  problem  of  dynamics  of  a  projectile  In  a  concentric  flexible  moving  tube. 
The  motion  of  the  tube  Is  accounted  for  by  using  the  Klrchhoff-Clebsch  theory 
of  deformed  rods.  However,  the  effect  of  projectile  balloting  Is  Ignored  and 
no  numerical  results  are  shown. 

This  report  presents  a  dynamic  analysis  of  a  balloting  projectile  In  a 
moving  gun  tube.  A  six  degree  of  freedom  mathematical  model  for  the 
projectile/gun-tube  system  Is  resorted  to  and  the  nonlinear  differential 
equations  of  motion  derived  using  a  Lagrangian  approach.  The  effects  of 
obturator  flexibility,  and  projectile  Impact  with  the  bore  at  the  bourrelet  or 
at  the  obturator  along  with  subsequent  rebound  are  included.  The  nonlinear 
equations  of  motion  are  solved  using  Gaussian  elimination  and  the  constant- 
average-acceleration  Integration  scheme.  Finally,  numerical  results  for  some 
sample  test  cases  are  obtained  and  discussed. 


5S.H.  Chu  and  F.K.  Soeohting ,  "Transverse  Motion  of  an  Accelerating  Shell," 
Picatinny  Arsenal ,  Dover,  NJ,  PA-TR-4314,  1972  \AD§  894S72). 

6 

F.J.  Ferdreaville,  "Analysis  of  the  Lateral  Motion  of  an  Unbalanced 
Projectile  in  a  Rigid  Gun  Tube,"  Sandia  Laboratories,  Albuquerque,  NM,  87115, 
SARD  74-0361,  1974. 

7 

F.J.  Perdreaville ,  "Arialysis  of  the  Lateral  Motion  of  an  Unbalanced 
Projectile  in  an  Elastic  Gun  Tube,"  Sandia  Laboratories,  Albuquerque,  NM, 
87115,  SAND  74-0362,  1974. 

8 

R.L.  Langhaar  and  A.P.  Boresi,  "Dynamics  of  a  Projectile  in  a  Concentric 
Flexible  Tube, "  BLM  Applied  Mechanics  Consultants ,  Contract  Report  No. 
ARBRL-CR-00501 ,  February  1983. 
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III.  SYSTEM  DESCRIPTION 


The  mathematical  model  employed  Is  a  slx-degree-of-freedom  coupled  model 
of  the  gun-tube  and  projectile  system.  An  Inertial  reference  frame  XYZ  Is 
selected  with  the  gun-tube  center  of  mass  at  (Xg,  Yg).  The  bore  may  be 
displaced  from  Its  original  position  due  to  recoil  and  vibration  with  the 
displaced  position  described  by  (Xg,  Yg,  e)  as  shown  In  Figure  1.  A  typical 
projectile  Is  represented  In  Figure  2.  The  origin  of  the  body-fixed  axes 
(x,y,z)  Is  fixed  on  the  bore  centerline  as  shown  In  Figure  3.  The  motion  of 
the  projectile  Is  described  by  the  coordinates  (Xp,  yp)  of  Its  center  of  mass 
relative  to  the  gun  and  the  yawing  angle  a  about  the  orthogonal  axis  2.  The 
projectile  Is  assumed  to  have  no  center  of  mass  offset  and  to  be  in  dynamic 
balance  about  Its  principal  axes.  Because  of  recoil  and  vibration,  the  gun 
bore  can  have  axial  as  well  as  lateral  and  rotational  motion.  As  the  pro¬ 
jectile  moves  down  the  gun  bore,  the  bourrelet  as  well  as  the  obturator  can 
impact  with  the  gun  bore  and  then  rebound  from  it.  This  effect  can  be  Includ¬ 
ed  by  using  simple  restoring  springs  to  represent  the  Interaction  between  the 
bore  and  the  projectile.  Both  the  gun  tube  and  projectile  are  assumed  to 
behave  as  rigid  bodies  and  rifling  of  the  projectile  Is  not  considered. 


IV.  THE  SYSTEM  KINETIC  ENERGY 

The  total  kinetic  energy  of  the  projectlle/gun-tube  system  Is 
T  *  Tgun  +  ^projectile 

=  *  m9  (S9?  *  %2)  *  *  !g 

*  *  mp  <*p2  *  fp2)  *  *  'p  (1) 

where  nig,  mp  represent  the  masses  of  the  gun  tube  and  projectile  and  Ig,  I 
represent  the  mass  moments  of  Inertia  about  axes  orthogonal  to  the  plane  of 
motion  at  the  respective  centers  of  gravity.  Expressions  for  the  quantities 
ftp  an<*  are  ^er^ve^  1°  Appendix  A. 
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Inertial  Frame  of  Reference 


Figure  1.  Displaced  Position  of  Gun  Bore 


Figure  2.  Projectile  with  Bourrelet  and  Obturator  within  Gun  Bore 


Figure  3.  Body-Fixed  Axes  and  Projectile  Displacements 


V.  BOURRELET/ GUN-BORE  CONTACT  HODEL 

When  the  projectile  contacts  the  gun  bore,  the  rebound  force  at  the 
bourrelet  can  be  Included  with  the  help  of  a  simple  spring  deflection  model, 
as  shown  In  Figure  4,  representing  the  stiffness  characteristics  of  the 
metallic  part  of  the  bourrelet  which  comes  Into  contact  with  the  gun  tube  upon 
Impact.  In  this  situation,  the  contribution  to  the  potential  energy  of  the 
system  would  be 

vbc  ■  *  kbc  4  <2> 

where  kbc  Is  the  appropriate  spring  constant  and  «bc  represents  the  displace¬ 
ment  of  the  projectile  Into  the  gun  bore,  that  Is,  the  distance  from  the  bore 
wall  to  the  undistorted  location  of  the  bourrelet  as  determined  by  the  orien¬ 
tation  of  the  projectile  axis.  This  displacement  Is  (see  Appendix  B) 

{bc  ■  |*P  +  lp  s,n  “I  -  Re1 .be  <3> 

where  represents  the  distance  between  the  projectile  CG  and  the  plane  of 
the  bourrelet  and  Rcl  br  Is  the  radial  clearance  between  the  bourrelet  and  the 
bore.  Contact  between  the  bourrelet  and  the  bore  occurs  when 


Figure  4.  Gbturator  and  Bourrelet/Gun-Bore  Contact  Models 
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This  contact  also  contributes  a  friction  force  applied  to  the  projectile  In 
the  negative  x  direction,  Its  magnitude  being 

Ffbc  ’  yfac  kbc  6fcc 

where  Is  the  coefficient  of  friction  for  the  bourrelet/gun-bore  contact. 

VI.  OBTURATOR  MODEL 


An  obturator  model  can  be  built  Into  the  balloting  analysis  by  including 
transverse  springs  as  shown  In  Figure  4,  In  the  plane  of  the  obturator, 
simulating  the  properties  Inherent  to  It.  The  spring  k  represents  the 

i  0 

stiffness  of  the  plastic  band  while  kQ  represents  the  stiffness  of  the  metal¬ 
lic  part  of  the  obturator  which  should  be  Included  In  the  obturator  stiffness 
when  the  metal  Impacts  the  gun  bore.  The  restoring  forces  in  these  springs 
will  tend  to  bring  back  the  projectile  CG  towards  the  centerline  of  the  gun 
bore  In  Its  motion  down  the  barrel.  The  spring  stiffness  values  are  to  be 
obtained  either  experimentally  or  analytically.  The  potential  energy  contri¬ 
bution  due  to  these  transverse  springs  will  be 


vobt  ■  *  ko  so  <"°  1n,pact> 

(6) 

V0bt  -  1  ko  so  *  1  ko  (5o  '  "cl  /  (upon  1mpact) 

I 

where  6Q  Is  the  projectile  displacement  at  the  obturator,  kQ  and  kQ  represent 
the  stiffness  of  the  plastic  band  and  metallic  part  respectively  and  Rr1  is 
the  radial  clearance  between  the  obturator  and  the  bore.  The  displacement  6Q 
Is  given  by 


yp  -  iQ  sin  a  (7) 

where  iQ  Is  the  distance  between  the  projectile  CG  and  the  plane  of  the 
obturator.  Contact  between  the  metallic  part  of  the  obturator  and  the  gun 
bore  occurs  when 


(8) 
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VII.  THE  SYSTEM  POTENTIAL  ENERGY 


The  total  potential  energy  of  the  system  Is  given  by 


V  -V+V  +  V  +V+V 

System  Gun  Proj  Bourrelet  Contact  Obt  Found  Moment 


■  9  Yg  +  fflp  9  Yp 

.  _  1  n 
+  I  -  a_  a 
_  n  n 
n 


*  *  kbc  4bc  +  *  ko  so  +  *ko  (4o  '  Rd.o>‘ 


(9) 


where  the  third  term  Is  to  be  Included  only  when  there  Is  contact  between  the 
bourrelet  and  the  bore,  the  fifth  term  Is  to  be  Included  only  when  the  metal¬ 
lic  part  of  the  obturator  Impacts  the  gun  bore,  and  the  last  term  Is  a  power 
series  due  to  the  foundation  moment  which  Is  a  large  nonlinear  effect  generat¬ 
ed  by  the  projectile  as  It  wobbles  down  the  gun  bore  and  which  can  be  measured 
experimentally.  The  quantities  mg  and  mp  represent  the  gun  mass  and  the 
projectile  mass  respectively  and  n  denotes  the  number  of  terms  In  the  power 
series  representation  of  the  foundation  moment,  which  Is  the  resistance  to 
wobbling  provided  by  the  stiffness  of  the  obturator  drive  band. 


VIII.  EQUATIONS  OF  MOTION 


Q 

Lagrange's  equations  of  motion  are 

|f  (4~)  “  4-  "  Qr  .  (r  -  1,  Z,  ...6)  (10) 

aqr  qr  nc 

where  L  Is  the  Lagranglan  function,  Q  Is  the  non-conservative  generalized 

rnc 

force  acting  In  each  coordinate  direction  qf  and  qf  denotes  the  generalized 
velocity  associated  with  the  generalized  coordinate  qp.  The  Lagranglan  L  Is 
the  difference  between  the  system  kinetic  and  potential  energies,  and  Is,  In 
this  case 


9 

L.  Meirovitah ,  "Analytical  Methods  in  Vibrations, "  Macmillan,  New  York,  1969 , 
pp.  30-50. 
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(11) 


l  ■  T-V  ■  {»  Bg  (I2  ♦  t2)  ♦  i  Ig  e2 

’  *  Ip ;2> 

•  «»,  9  Ig  *  rnp  9  vp  ♦  J  kbe  a2p  *lk0«I 

♦  *  k'  (J„  -  R,,  ,)!  *  I  J  a.  a") 

oo  ci »o  ^  n  n 
n 


where  the  quantities  flp  and  f  are  derived  In  Appendix  A.  dbc  is  given  In 
Appendix  B  and  iQ  and  Q  are  as  defined  In  the  section  on  Obturator  Model. 

Application  of  Lagrange's  equations  yields  the  following  nonlinear, 
coupled  equations  of  motion  for  the  projectile/ gun-bore  system 

•4  44  44 

ql  ■  xg:  (mg  +  V  X9  "  mp  *xp  s1n  9  *  yP  C0S  6>  9  +  mp  xp  cos  6 

44  O  O 

-  rn  y„  sin  6  *  mj^  cos  9  +  2m„x„0  sin  0  -  m  v  9^  sin  9 
P  P  P  P  P  P  FP 

+  2  ">^0  cos  9  (12) 

44  44  44 

q,  *  Y  :  (m„  +  mj  Y„  +  in  (x„  cos  0  -  y„  sin  9)  0  +  m„  x„  sin  0 

c  g  g  p  g  p  p  'p  p  p 

+  "L  y„  cos  0  ■  mj<  02  sin  e  -  2mJLe  cos  0  +  m^yj2  cos  0 
P  P  P  P  P  P  rP 

+  2  sin  0  -  (mp  +  mg)g  (13) 


^3  -  0:  mp(xp  sin  e  +  yp  cos  0)  Xg  +  mp  (xp  cos  0  -  yp  sin  e)  Yg 

o  O’*  ••  •• 

*  (Ig  ♦>p'p*'pi,)*"p  vp  Xp  *  «„  xp  yp 

*  ”p  (2  5  <xpxp  *  yp  Vp)  -  *5  Clip  s'"  «  +  VP  COS  e) 

+  0  (xp  cos  0  -  yp  sin  0)3  ♦  fg  C(x  cos  0  -  yp  sin  0) 

-  0  ( xp  sin  0  +  yp  cos  0)]}  -  mp  ([*g  yp  -  fg  xp)  0 

-  <typ  +  Vp)]  s,n  6  '  [(Vp  *  Vp1  5 

+  *  M„)]  C0S  0>  ’  0 

g  p  g  p 


(14) 
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«l4  "  xp:  mp  <Xg  +  *g  *  •  [Xp(cos  e  -  sin  6)  -  yp(s1n  0  +  COS  0)3 

+  (sin  0  +  cos  0)  Xp  +  (cos  e  -  sin  0)  yp  -  2  yp§  (sin  0  +  cos  0) 

A  A 

+  2  xp0  (cos  0  -  sin  0)  -  (Xp+  yp)  h  cos  0  *  (yp  -  xp)  0  sin  0 

-  Xp  0^  -  yp  0  -  (fgCOS  0  -  itg  sin  0)  0  +  g  sin  0} 

-  fx(t)  -  ubc  kbc(|yp  ♦  ip  «|  -  Rc1§bc)  sgn  (xp)  (15) 

q5  "  yp:  mp  {"  *gs  fn  8  +  Yg  cos  0  +  xp  0  +  yp  ♦  2  xp  0  -  tg  0  sin  0 

-  X  5  cos  e  -  y  S2  *  5  (J  cos  0  +  t  sin  0) 

9  r  9  9 

*  9  «s  »>  *  kbc  <V  *P  “  '  Pcl  ,bc  s9n  (yp  ♦  ‘p  «>> 

*  <ko  *  ki)  (yp  •  V>  •  Rcl,o  ko  S3"  (yp  -  l„  a)  ■  fy(t)  (I6) 

"6  ’  a:  'p  •  *  kbc  (yp  *  *pa  *  "cl  .be  s9n  (yp  *  ‘p  «»  ‘p 
*  (k0+  kg)  t0  (yp  -  t0  a) 

+  Rcl  0  ko  *o  sgn  *yp  “  V°  +  1  an  °nl  +  cl“  +  c2  a  U|-  m  (t)  (17) 
Equations  (12-17)  can  be  written  down  In  the  form 

[M]{q>  +  [C]  {q}  +  [K]  (q)  •  {F}  (18) 

where  [M],  [C],  and  [K]  are  6x6  matrices,  (F>  Is  a  force  column  matrix  and 
(q>  Is  a  vector  of  generalized  coordinates.  The  elements  of  these  matrices 
are  as  follows: 

Mn  "  (mg  +  V  •  m13  ■  *  mp  (*p  s1n  0  +  yp  cos  6)  , 

M14  *  mp  cos  e  *  M15  "  ‘  mp  sin  0  , 

M22  »  (mg  +  mp)  ,  M23  -  mp  (xp  cos  e  -  yp  sin  0)  , 

M24  "  mp  s*n  6  ’  M25  *  mp  cos  0  » 

M31  *  '  fflp  (xp  s1n  9  +  yp  cos  * 

M32  *  mp  *xp  cos  9  "  yp  s1n  9^  * 


17 


(19) 


*33  *  ’g  *  Vxp  +  yp’  •  *34  '  *  Vp 

*35  "  mpxp'  •  *41  *  mp  •  *42  ’  mp  • 

*43  mp  [x  (cos  e  -  sin  e)  -  yp  (sin  8  +  cos  8)]  » 

M44  *  mp  6  +  cos  e)  •  M45  *  mp  (cos  6  “  s^n  0)  * 

Mgi  ■  -  mp  sin  8  ,  Mgg  »  mp  cos  e  , 

M53  “  mp  *p  •  M55  "  mp  ’  M66  "  *p  * 

»  -2  mp  8  sin  8  ,  Cjg  ■  -2  mp  8  cos  6  , 

C23  “  2  mp  (x  cos  8  -  yp  sin  8)  , 

C31  *  -  {(ip  sin  8  +  yp  cos  8)  +  8  (x  cos  8  -  yp  sin  0)>  mp  , 

C32  ■  {(xp  cos  e  *  yp  sin  8)  -  0  (xp  sin  0  ♦  yp  cos  8)}  mp  , 

c33  ■  -  »p  {(*g  3P  -  fyp>  S'"  6  -  <»g*p  +  Vp5  C”  6>  ’ 

C34  ■  "Ip  (*g  s'"  «  -  fg  COS  9)  . 

^35  “  mp  *^g  s1n  0  +  s*n  0^  * 

C41  *  mp  ®  s1n  0  •  c42  *  -  mp  ®  cos  8  ,  (20) 

C43  c  -  mp  (xp8  ♦  yp)  ,  C44  -  2mpe  (cos  8  -  sin  8)  , 

C45  *  *20  mp(s-|n  8  +  cos  8)  , 

Cgj  ■  -mp  e  cos  8  ,  Cg2  ■  -«np  0  sin  0  , 

C53  “  mp  (*g  cos  0  +  *g  sin  e>  *  c54  "  2V  * 

C66  "  C1  +  c2  !iL 
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(21) 


*14  *  “mpSZ  cos  e  •  *15  "  mp®2  s1n  6  • 

K24  "  'V2  s1n  0  »  *25  *  *mp®2  cos  •  . 

K34  *  ^p  *  *p  •  *35  “  V  ip  • 

*44  "  ""V®2  ^C0S  ®  +  s1n  ®^  *  *45  "  "’ll®2  ®  *  C0S  ®^  * 

K55  *  'V2  +  kbc  *  ko  +  K  • 

he  ■  "be  ‘p  -  <ko  *  ki>  *0  • 

k65  ■  “be  ‘p  -  <k0  *  ki>  ‘o  - 

k66  *  kbc  ‘p  *  <ko  *  ki>  ‘o  *  S  *n  a"‘?  • 

F2  ■  -(mp  +  mg)  g  , 

f4  *  ""p  3  s,n  '  *  fx(t)  -  “be  kbc  (k  *  lp“l  *  "cube1  s«n  <V  ’ 

1  (22) 
F5  ■  -mp  9  cos  •  ♦  kbc  Rc,ibc  sgn  (yp  ♦  *p«>  *  fy(t) 

*  "cl.o  K  S9n  %  '  V>  ’ 

f6  ■  mCO  *  k„c  Pc1<bc  »„  sgn  (y„  ♦  y)  •  »clt„  k’  »„  sgn  (yp  -  «„»)  ■ 
All  other  elements  are  zero. 


IX.  SOLUTION  TECHNIQUE 

The  equations  of  motion  (Equations  18)  for  the  projectll e/gun  tube  system 
are  coupled,  nonlinear  differential  equations  which  can  be  solved  only  numer- 

e 

Ically.  If  the  values  of  q's  and  q's  are  known  at  any  Instant,  their  substi¬ 
tution  In  [M],  [C],  [K]  and  will  lead  to  a  set  of  6  simultaneous  algebraic 
equations  with  the  q's  as  the  unknowns.  Elements  of  these  matrices  are  first 
computed  from  given  Initial  values  of  q's  and  q's  at  the  first  time  step  and 
then  the  simultaneous  equations  generated  are  solved  by  Gaussian  elimination 
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to  find  the  q's  which  lead  to  the  solution  for  the  subsequent  time  step. 
Essentially,  we  have  a  set  of  six  coupled  nonlinear  differential  equations  to 
be  solved  at  ekch  time  step.  The  solution  to  these  equations  may  be  obtained 
by  Newmark's  constant  average-acceleration  method  of  Integration10  which  Is  a 
self-starting  technique  that  yields  a  step-by-step  solution  up  to  a  specified 
time  limit,  once  the  Initial  positions  and  velocities  are  supplied.  The  time 
step  can  be  monitored  during  the  process  of  computation  If  needed.  At  Impact 
end  rebound  between  the  projectile  and  the  gun  bore,  either  at  the  bourrelet 
or  at  the  obturator,  the  solution  of  the  equations  of  motion  need  not  be 
Interrupted  although  additional  forces  will  come  Into  play  from  the  activation 
of  springs  kbc  and  k^  and  from  the  frictional  effect  between  the  bourrelet  and 
the  bore  when  the  bourrelet  slides  on  the  walls  of  the  gun  tube.  The  projec- 
tlle/gun  tube  dynamic  system  response  can  then  be  studied  at  any  desired  time 
and  time  histories  of  all  coordinates  describing  the  system  behavior  can 
easily  be  obtained. 

The  first  consideration  In  selecting  a  numerical  Integration  scheme  for 
multi  degree  of  freedom  systems,  as  Craig11  points  out,  should  be  Its  stabil¬ 
ity.  In  most  cases.  It  Is  desirable  to  use  a  method  that  Is  unconditionally 
stable  such  as  the  constant-average-acceleratlon  method.  This  method  also 
produces  no  amplitude  error,  that  Is,  there  Is  no  numerical  dissipation, 
regardless  of  the  time  step  size.  The  period  error,  too.  Is  small.  However, 
for  some  multi  degree  of  freedom  systems,  It  is  desirable  to  have  numerical 
dissipation  to  filter  out  the  response  of  less  accurate  higher  modes,  which  Is 
similar  to  truncating  higher  modes  In  a  mode-superposition  solution.  Accord¬ 
ing  to  Craig,  Newmark's  linear  acceleration  method  has  an  advantage  over  his 
constant-average-acceleration  method  In  that  It  provides  this  damping  for 
higher  modes.  One  difficulty,  however.  Is  that  it  Is  not  unconditionally 
stable.  In  order  to  ensure  an  accurate  solution  to  the  balloting  projectile 
problem,  both  of  these  methods  were  resorted  to  and,  for  small  time  steps,  no 


Bathe  and  E.L.  Wilson ,  "Numerical  Methods  in  Finite  Element  Analysis," 
Prentice  Hall ,  Englewood  Cliffs,  N.J.,  1976,  pp.  322-326. 

Hr.R.  Craig,  Jr.,  " Structural  Dynamics  -  An  Introduction  to  Computer 
Methods,"  John  Wiley,  New  Fork,  1981,  pp.  461-463. 
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difference  In  the  results  was  noted.  -  It  was  therefore  decided  to  use  the 
constant-average-acceleratlon  method  for  this  work  because  of  Its  uncondition 
al  stability. 


X.  COMPUTER  CODE  DESCRIPTION 

All  computer  code  used  In  this  study  was  written  In  Pascal  on  the 
Hewlett-Packard  9el6  Series  200  Desk  Top  Computer.  The  main  program,  NEWMARK, 
performs  numerical  Integration  based  on  either  the  linear  acceleration  or 
constant-average-acceleratlon  method.  MOPBOX  (Matrix  Operations  tool  BOX) 
provides  various  low-level  routines  to  manipulate  matrix  and  vector  entitles. 
Including  addition,  subtraction,  multiplication,  etc.  There  Is  also  a  routine 
that  solves  linear  algebraic  equations  by  Gauss-Jordan  Elimination.  Graphics 
support  Is  provided  by  the  PLOTPAC  Graphic  Plotting  Package.  For  a  complete 
listing  of  the  code,  see  Appendix  D  -  Computer  Program  Listing. 

NEWMARK,  as  well  as  Its  supporting  modules,  MOPBOX  and  PLOTPAC,  were 
tested  using  the  example  of  a  linear,  two  degree  of  freedom  system  given  In 
Bathe  4  Wilson10  on  pp.  324-325.  The  Input  used  Is  shown  on  page  76,  Appendix 
D.  Results  from  NEWMARK  were  In  exact  agreement  with  those  given  In  this 
reference. 


XI.  NUMERICAL  RESULTS  AND  DISCUSSION 

The  balloting  projectile  model  and  analysis  developed  have  been  utilized 
to  simulate  the  dynamical  behavior  of  an  APFSDS-T  (Armor  Piercing  Fin  Stabil¬ 
ized  Discarding  Sabot  with  a  Tracer)  projectile  In  a  120  mm  M1A1  Tank  Main  Gun 
Tube  In  the  following  situations: 

Case  1  Projectile  CG  at  the  Obturator  (iQ  »  0) 

(a)  Initial  condition:  yQbt  ■  0  ,  a  -  0.054  deg  ,  correspond¬ 
ing  to  a  maximum  tilt  of  the  projectile  with  no  displacement  at 
the  obturator. 
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(b)  Initial  condition:  yQbt  •  -.0026  In  ,  a  ■  0.054  deg  , 
corresponding  to  same  projectile  tilt  as  In  (a)  In  conjunction 
with  a  displacement  at  the  obturator  equal  to  one-half  of  the 
radial  clearance  between  the  projectile  and  the  gun  bore  at  the 
obturator. 


Case  Z  Projectile  CG  Between  the  Obturator  and  the  Bourrelet.  2  Inches 

from  Obturator  [iQ  ■  c  incnes) 

(a)  Same  as  In  1(a). 

(b)  Same  as  In  1(b). 

Other  data  used  Is  as  follows: 

Gun  weight  ■  4000  lbs 

Projectile  weight  ■  16  lbs 

Projectile  cross-sectional  area  •  17.53  In 


■  0.583  1b-1n-sec8 

■  150.0  lb-1n-sec8 

■  71  lb-ln-sec,  approximately  5%  of 
critical 


-  0 

*  5.0  x  107  lbs/in 
-  5.0  x  106  lbs/in 
■  1.0  x  108  lbs/in 


Rcl,  be 
Rcl.  o 


■  5.5  Inches 
B  0.0052  Inches 
*  0.0052  Inches 


Time  history  of  gas  pressure  p(t)  on  the  projectile  as  given  In  Figure  5. 
The  curve  for  p(t)  shown  Is  a  modification  that  results  from  Inclusion  of 
the  effect  due  to  friction  at  the  obturator/bore  Interface. 


Projectile  foundation  moment  as  given  In  Figure  6. 
The  time  step  used  for  all  solutions  Is  0.005  msec. 
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Gun  Horizontal  and  Vertical  Displacements  and  Projectile  x-Dlsplacement 
Response 

The  time  histories  of  the  horizontal  gun  displacement  Xg(t),  the  vertical 
gun  displacement  Yg(t)  and  the  x-dlsplacement  xp(t)  of  the  projectile  generat¬ 
ed  by  the  analysis  are  shown  In  Figures  7-9.  These  are  similar  for  all  the 
cases  analyzed  and  no  significant  changes  are  noted  from  one  case  to  the  next. 


Figure  7.  Time  History  of  Horizontal  Gun  Displacement  (Xg) 


Time  (msec) 

Figure  8.  Time  History  of  Vertical  Gun  Displacement  ( Y  ) 
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350 


Figure  9.  Time  History  of  X-Dlsplacement  of  Projectile  (Xp) 

Gun  Angular  Displacement  Response 

The  time  histories  of  the  angular  displacement  © ( t)  of  the  gun  for  the 
various  cases  are  represented  In  Figures  10-13.  For  a  projectile  with  Its  CG 
located  at  the  obturator,  and  Initially  cocked  In  the  bore  with  no  displace¬ 
ment  at  the  obturator,  the  angular  gun  displacement  e(t)  gradually  Increases 
with  time  as  shown  In  Figure  10  until  the  projectile  exits  the  bore  with  no 
Impact  occurring  between  the  bourrelet  and  the  bore  at  any  time  at  all.  When 
an  Initial  vertical  displacement  Is  Imposed  at  the  obturator,  oscillations  of 
the  tube  result  from  continual  Impact  between  the  bourrelet  and  the  bore  as 
shown  In  Figure  11.  For  the  case  In  which  the  projectile  CG  Is  located  away 
from  the  obturator  and  the  projectile  is  Initially  cocked  In  the  bore  with  no 
displacement  at  the  obturator,  the  gun  angle  increases  until  the  first  Impact 
between  the  bourrelet  and  the  bore  and  then  begins  to  oscillate  as  seen  In 
Figure  12.  With  an  additional  Initial  vertical  displacement  Imposed  at  the 
obturator,  the  oscillations  are  seen  to  be  small  until  the  first  moment  of 
Impact  between  the  bourrelet  and  the  bore,  following  which  oscillation  ampli¬ 
tudes  Increase  with  further  Impacts. 


0.0000125 


Figure  10.  Time  History  of  Gun  Angular  Displacement  (e)  for  Case  la 
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Figure  11.  Time  History  of  Gun  Angular  Displacement  (e)  for  Case  lb 


Projectile  y-Dlsplacement  Response 

For  case  1(a),  Figure  14  essentially  reveals  a  low-frequency  type  oscil¬ 
latory  behavior  that  results  from  the  projectile  not  Impacting  the  gjn  bore  at 
all.  For  case  1(b),  owing  to  continual  Impact  between  the  bourrelet  and  the 
bore,  the  oscillations  are  seen  to  be  of  higher  frequency  as  shown  In  Figure 
15.  For  case  2(a),  Figure  16  shows  a  gradual  fall  of  displacement  until  the 
occurrence  of  the  first  Impact  between  the  bourrelet  and  the  bore  leading  to 
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Figure  14.  Time  History  of  Y -Displacement  of  Projectile  (yp)  for  Case  la 
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Figure  15.  Time  History  of  Y-Olsplacement  of  Projectile  (yp)  for  Case  lb 


Projectile  Yaw  Response 


The  time  histories  of  the  yawing  motion  of  the  projectile  a(t)  generated 
for  the  various  cases  are  represented  In  Figures  18-21 .  For  case  1(a),  the 
response  Is  a  low  frequency,  sinusoidal  type  of  motion  as  shown  In  Figure  18, 
which  results  from  the  projectile  not  Impacting  the  gun  tube  at  all.  For  case 
1(b)  (see  Figure  19),  because  of  continual  Impact  between  the  bourrelet  and 
the  bore,  oscillations  of  larger  frequency  occur.  For  case  2(a)  (see  Figure 
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Figure  16.  Time  History  of  Y-D Is placement  of  Projectile  (yp)  for  Case  2a 
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Figure  17.  Time  History  of  Y-Dlsplacement  of  Projectile  (yp)  for  Case  2b 


20),  the  bourrelet  Is  initially  In  contact  at  the  top  of  the  bore.  As  the 
projectile  travels  down  the  gun  barrel,  the  bourrelet  drops  making  Impact 
twice  before  exiting.  For  case  2(b)  (see  Figure  21),  pre-impact  oscillations 
are  relatively  small,  but  become  large  following  bourrelet  contact. 
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Figure  19.  Time  History  of  Yawing  Motion  of  Projectile  (a)  for  Case  lb 
Effect  of  Using  a  Smaller  Time  Step  for  Integration 

In  order  to  determine  the  effect  of  a  smaller  time  step  on  the  dynamic 
behavior  of  the  system,  case  2a  Is  also  analyzed  using  At  »  .0005  msec.  There 
Is  no  appreciable  change  In  results  except  that  the  oscillation  period  becomes 
slightly  smaller,  which  Is  In  agreement  with  Craig's  observation1*. 
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Figure  20.  Time  History  of  Yawing  Motion  of  Projectile  (a)  for  Case  2a 
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Figure  21.  Time  History  of  Yawing  Motion  of  Projectile  (a)  for  Case  2b 


Other  Observations 

The  results  shown  Indicate  that  the  projectile  motion  oscillates  about 
the  limits  set  by  the  clearance  between  the  projectile  bourrelet  and  the  gun 
bore.  The  maximum  longitudinal  deceleration  of  the  gun  is  220  g's  and  occurs 
at  about  2  msec  for  all  the  cases  analyzed.  The  maximum  lateral  acceleration 
or  deceleration  of  the  gun  Is  lowest  (1  g)  for  Case  la,  in  which  the  projectile 
CG  Is  at  the  obturator  and  the  projectile  Is  Initially  cocked  In  the  bore  with 
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no  vortical  displacement  at  the  obturator.  It  Is  largest  (14  g)  for  Case  2b, 
In  which  the  projectile  CG  Is  2  Inches  away  1?rom  the  obturator  and  the  pro¬ 
jectile  Is  Initially  cocked  In  the  bore  with  a  vertical  displacement  at  the 
obturator.  The  maximum  projectile  longitudinal  acceleration  is  seen  to  be 
57,000  g,  the  same  for  all  cases.  The  maximum  projectile  lateral  acceleration 
or  deceleration  Is  the  lowest  (45  g)  for  Case  la  and  highest  (48,000  g)  for 
Case  2b.  The  maximum  projectile  yaw  rate  Is  the  lowest  (0.8  rads/sec)  for 
Case  2a  and  highest  (9.5  rads/sec)  for  Case  2b,  with  the  maximum  yaw  being  of 
the  order  of  .001  radians  for  all  cases.  Details  are  shown  In  Table  1.  In 
general,  the  projectile  whose  CG  Is  situated  away  from  the  obturator  and  which 
Is  Initially  cocked  In  the  bore  with  a  prescribed  vertical  displacement  at  the 
obturator  (Case  2b)  seems  to  exhibit  the  most  vigorous  dynamic  behavior,  as 
expected. 

For  the  projectile  whose  CG  Is  at  the  obturator  and  which  Is  Initially 
cocked  In  the  bore  but  has  no  prescribed  Initial  vertical  displacement  at  the 
obturator,  the  results  Indicate  that  there  Is  no  Impact  at  all  between  the 
bourrelet  and  the  bore  throughout  the  In-bore  period.  Consequently,  Ecuatlon 
(17)  which  Is  the  equation  of  motion  corresponding  to  the  sixth  generalized 
coordinate,  namely,  the  projectile  yawing  motion,  reduces  to  an  uncoupled 
second  order  differential  equation  In  the  coordinate  a.  If  the  foundation 
moment  curve  Is  linear  over  the  range  of  the  In-bore  period,  as  is  the  case 
here,  this  equation  of  motion  further  reduces  to  a  linear  uncoupled  second 
order  differential  equation  In  o.  Thus,  so  long  as  the  projectile  CG  Is 
situated  at  the  obturator  and  the  obturator  has  no  Initial  vertical  displace¬ 
ment,  the  projectile  could  be  assumed  to  be  a  single  degree  of  freedom  system 
for  obtaining  a  good  Initial  estimate  of  yawing  response.  However,  for  other 
situations,  a  more  rigorous  multi  degree  of  freedom  system  analysis  would  seem 
necessary,  as  the  results  of  this  work  Indicate. 


XII.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 

In  order  to  ensure  accuracy  In  Internal  ballistics,  a  realistic  predic¬ 
tion  of  the  yawing  motion  of  a  projectile  Is  necessary.  A  dynamic  analysis 
incorporating  a  six  degree  of  freedom  mathematical  model  of  a  projectile/gun- 
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TABLE  1.  COMPARISON  OF  MAXIMUM  RESPONSES  FOR  THE  DIFFER^  CASES  ANkLYZEO 
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•-tube  system  Is  presented  In  this  report.  The  nonlinear  differential  equations 
of  motion  of  the  system  are  derived  using  a  Lagranglan  formulation  with  the 
effects  of  obturator  flexibility  and  projectile  Impact  and  rebound  at  the 
bourrelet  and  the  obturator  Included.  Using  Gaussian  elimination  and 
;  ‘4Hewnark 1 s  constant-average-acceleratlon  Integration  scheme,  the  nonlinear 
equations  are  solved  to  yield  usable  dynamical  response  information  for  some 
specific  test  cases. 

The  analysis  and  results  Indicate  that  as  long  as  the  CG  of  the  projec¬ 
tile  Is  located  at  the  obturator  and  no  Initial  vertical  displacement  Is 
prescribed  at  the  obturator,  the  projectile  can  be  treated  as  a  simple  single 
degree  of  freedom  system  for  estimating  Its  yawing  response.  For  more  complex 
situations,  however,  a  more  rigorous  analysis  employing  a  larger  number  of 
degrees  of  freedom  such  as  the  one  suggested  In  this  report  seems  to  be  In 
order. 

The  following  recommendations  for  further  work  are  bound  to  constitute 
additional  contributions  to  the  general  problem  of  the  dynamics  of  a  balloting 
projectile  In  a  moving  gun  tube. 

1.  Account  for  the  three-dimensional  motion  of  the  projectile  and  the 
gun  employing  Euler's  angles. 

2.  Treat  the  projectile  as  a  rigid  body  moving  within  an  elastic  gun 

tube. 

3.  Investigate  the  effect  of  considering  the  projectile  as  a  flexible 

body. 

4.  Include  the  effect  of  rifling  of  the  projectile. 
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APPENDIX  A 


PROJECTILE  CENTER  OF  GRAVITY  DISPLACEMENTS  AND  VELOCITIES 


The  position  vector  ?  of  the  projectile  center  of  mass  In  the  Inertial 
frame  of  reference  can  be  written  as 

%  "  *p/g  +  *G  (A’1} 

when  rp,6  represents  the  position  vector  of  the  projectile  center  of  mass  with 
respect  to  the  gun  tube  center  of  mass(  and  fQ  represents  the  position  vector 
of  the  gun  tube  center  of  mass  In  the  Inertial  frame  of  reference.  Using  1 

A  A  A 

and  j  as  unit  vectors  along  X  and  Y  axes  and  t  and  m  as  those  along  the 
body-fixed  x  and  y  axes,  we  can  write 


*P/G  *  V  +  *p  ■ 

.  A  A 

rr  s  Xn1  +  YJ 
»  9  9 

A  A  A 

where  i  ■  cos  e  1  +  sin  e  j 
m  *  -sin  e  1  +  cos  e  J 

Substituting  equations  (A-2)  Into  (A-l),  we  get 


(A-2) 


(A-3) 


■  Xp'  *  V  (A-4) 

where  X  -  (x  cos  e  -  y  sin  0  +  X  ) 

(A-5) 

Yp  •  (*p  *1"  0  +  yes  •  +  Yg) 

Equations  (A-5)  can  be  differentiated  to  yield  the  relative  velocities  Xp 
and  ?  for  the  projectile  CG  as  follows. 


Xn  •  -x  sin  e  0  +  x„  cos  e  -  y„  cos  e  0  -  y^sln  0  +  X 
P  P  P  P  P 

■  x„  cos  e  0  ♦  x.  sin  0  -  y„  sin  e  e  +  y.cos  e  + 

P  P  P  P  P  9 


(A-6) 
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APPENDIX  B 


PROJECTILE  DISPLACEMENT  AT  BOURRELET 

The  radial  displacement  vector  of  the  projectile  centerline  In  the  plane 
of  the  bourrelet,  measured  from  the  bore  centerline  Is  (see  Figure  22) 

r  ■  rB  *  r8* 

■  <Vp  +  Vs  *  V  -  <V/s +  V 

♦  ♦  4 

"  rB/P  rP/G  ’  rB*/G 


1  cos  «  1  +  sin  a  m  +  xA  +  m  -  t  t  (B-l) 

P  P  P  P  9 

A 

where  P  denotes  the  CG  of  the  projectile  and  G  denotes  the  CG  of  the  gun,  t 
and  m  are  unit  vectors  along  and  perpendicular  to  the  gun  axis,  and  1  and  j 
are  unit  vectors  along  the  Inertial  X  and  Y  axes.  Here  which  Is  the 
Instantaneous  distance  of  the  bourrelet  from  the  gun  CG  measured  along  the  gun 
axis,  can  be  seen  to  be 


(B-2) 

(B-3) 

(B-4) 


i  «  x  +  i  cos  o 
9  P  P 

Using  the  above  approximation,  It  can  be  seen  that 
r  •  A1  +  Bj 

where  A  ■  -(tp  sin  a  +  yp)  sin  e 

B  ■  (ip  sin  a  ♦  yp)  cos  9 
Contact  occurs  when 

H  5  "cl.bc  <8-5> 

where  R£^  Is  the  radial  clearance  between  the  bourrelet  and  the  gun  bore. 
The  displacement  of  the  projectile  Into  the  bore  Is  then 


’be 


cl  ,bc 


yn  +  sin  a 
P  P 


cl  ,bc 


(B-6) 
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Figure  22.  Bourrelet  Impacting  Bore 
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GENERALIZED  FORCES  ACTING  ON  SYSTEM 


The  virtual  power  developed  due  to  generalized  forces  through  Infinites- 
Imal  virtual  velocities  compatible  with  system  constraints  can  be  represented 
% 


ip 


tot 


6 

1  «k 


iv. 


k-1 


Qy  &X.  +  Qv  +  Q, 


V* 


■*■0  lx 

*P  P 


+  Qy  «yp  +  Qa 

p 


(C-l) 


where  Q.,  represents  the  nonconservative  generated  force  associated  with  the 

*  th 

generalized  velocity  vk  In  the  direction  of  the  k  generalized  coordinate. 

The  system  virtual  power  will  also  be  that  resulting  from  the  application  of 

given  nonconservative  forces  acting  on  the  system,  which  Is  the  sum  of  the 

following  contributions. 


6P  due  to  force  fx(t)  generated  by  the  gas  pressure  on  the  projectile 
•  fx(t)  S*p 

6P  due  to  force  fy(t)  on  projectile  •  fy(t)  iyp 

6P  due  tc  moment  m(t)  on  projectile  »  m(t)  ia 

6P  due  to  friction  at  contact  surface  between  the  bourrelet  and  the  bore 

’  Ffbc  ,5p  s,n  (V  '  *"bc  kbc  (|yP  *  v|  •  Rcl,bc>  %  S3n  <*p> 

for  small  o,  upon  contact. 

6P  due  to  quadratic  yaw  damping  »  -(c^a  +  a  |a| )  6a 


Thus-  ■  (fx(t)  -  pbc  kbc  (|yp  ♦  tp  «|  -«c1ibc)  59"  (*„»  5i, 


+  {fy (t)  6y  )  +  (m(t)  -  (Cj  a  +  c2  i  a  )} 


ia 


(C-2) 


A  comparison  of  equations  (C-l)  and  (C-2)  then  yields  the  following  gener¬ 
alized  forces 


*7 


(C-3) 


Qe  -0 
-  “be  kbc  < 


c,i  -  c2  i  |i| 

In  the  above,  m(t)  Is  any  external  yawing  moment  acting  on  the  projectile,  vbc 

Is  the  coefficient  of  friction  at  the  bourrelet  and  bore  contact,  c^  and  c2 

are  first  order  and  second  order  damping  coefficients  for  the  yawing  motion  of 

the  projectile,  and  f  (t)  and  f  ( t)  are  applied  forces  on  the  projectile  In  x 

x  y 

and  y  directions  given  by 

fx(t)  *  p(t)  Ap  COS  a  , 

(C-4) 

fy(t)  *  p(t)  Ap  Sin  a  , 

where  Ap  Is  the  area  of  the  projectile  cross  section  on  which  the  gas  pressure 
p(t)  acts. 


y_  +  A<* 
P  P 


-  W  S9n  ^p1 


\  ‘ 

9  9 

ix  ■  f,(‘) 

P 

Qyp  *  fy(t) 
Qa  •  m(t)  - 
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tucadl 

•sysprog* 

program  rtewmark  f Input  ,  output,  keyboard): 


<  . . . . . . .  ) 

(4  •  ) 

(*  N  E  W  H  A  R  K  *) 

(*  Numerical  Integration  Program  *> 

<•  *  ) 

<•  Author:  John  Baugh,  Bat  tel le-Northwest  Laboratories  •) 

(•  Data:  Novanber  1,  1985  *> 

<  •  •  > 

<*  Abstract:  This  program  determines  displacements,  velocities,  and  «) 

<  «  accalerat ions  using  *ha  Newmark  numerical  integration  method  for  •) 

<•  nonlinear,  multi-degree-of-freedom  systems.  Both  the  constant-  *) 

<•  average-acceleration  and  linear  acceleration  methods  are  available  •) 

<•  to  the  u#er.  The  system  matrices  <e.g.  mass,  stiffness,  etc.)  are  •  ) 

<•  defined  by  subroutines  in  the  file  'USER. TEXT’.  Subroutines  are  ’*  > 

<4  called  at  each  time  step  and  may  be  functions  of  the  current  •> 

<*  displacement,  velocity,  and  time  vectors.  •) 

(  *  *  > 

. . . . . . 444444 . . 


•search  ’M0P80X’#  <•  access  matrix  operations  •) 

import  mopboxi 

•search  'PLOTPAC’S  (•  access  plotting  routines  «) 

import  plotpac: 

const 

MAXOOFS  -  Gi 
PI  -  3. 141592654, 
type 

string80  «  s tr i ng[ 80 1 1 
plot_record  • 
record 

dlapl  :  plot_array: 
veloc  :  plot_arrayi 
accel  :  plot_array 
end: 

var 

outfile  :  texti 
filename  :  strtng80: 
ndofs  :  integer: 
npts  :  integer: 
ch  :  char: 
step_aize  :  real: 
maxtime  :  real: 
time  :  plot_array: 

data  :  array  II.. MAXOOFS]  of  plot_record: 

. .  USER-DEFINED  INPUT  FUNCTIONS  . . *  > 

•include  ’USER. TEXT' • 

. . .  UTILITIES  . . . . 

function  open_ou t f i le( f 1 lename  :  strino30>  :  boolean: 

(*  open  an  output  file  if  possible  •) 
begin 
try 

rewrl tel  out f i le ,  filename): 
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open_outfile  : «  truei 
recover 

ooen_outflle  ;■  false 
endi  ;«  open_out#ile  * ) 


function  ctoi  (ch  :  char)  :  intagen 

(•  convert  a  character  to  an  integer  digit  •) 

begin 

ctoi  :»  ord(ch)  -  ord( '0*  ) 
endi  (•  ctoi  *) 


procedure  output_deta  (npts  :  integer)! 

(*  output  the  displs,  vetoes,  and  accels  to  the  specified  file  •> 
var 

i  ,  dof  :  integer! 
begin 

wn  taint  out  f  i  le  ,  '  dof  tine  ( «S  >  displ 

' veloc  accel 1  > i 

wr i te 1 n( ou  t  f i le  )« 
for  i  1  to  npts  do 
begin 

for  dof  :■  1  to  ndofs  do 
begin 

ur i tel  out f i le ,  dof:4,  ’  ’>! 

wr  l  tel  out  f  i  le  ,  timet i J : 10:4  ,  *  Ml 
wr  i  te<  out  f  i  le  ,  dat  at  dof  1 .  di  spl  C 1  ] :  1 0: 4  ,  '  Ml 

wr  i  tel  out f i le ,  datat dof 1 . veloc t i  1 : 1 0: 4  ,  ’  *)i 

writelnloutfile,  do  tat  do f 1 .accelt i } : 1 0: 4  > 
end! 

wri telnl out  f 1 le ) 
end 

endi  (•  oulput_data  ♦) 


procedure  display_plots  vdof,  npts  :  integer)! 

<*  display  the  displ,  veloc,  and  acccl  plots  at  the  specified  dof  *) 
procedure  display  <*,  y  :  plot_arrayi  title  :  plot_str ing )  i 

vir 

ch  :  char; 
begin 

if  initgraphics  then 
begin 

plot!:*,  y  ,  rots,  1  ,  title/; 
readl K eyboard ,  ch); 
termgraphics 
end 

end ;  1  *  di splay  *  > 


pen :  n 

display. time  ,  datat cof  1 . displ  ,  'Displacement  vs  Time')i 
d: splayi t ime  ,  da t a [ do f 3 . ve 1 oc  ,  'Velocity  vs  Time'); 
di sp lay. t i me ,  datat aof 3 .accel ,  'Acceleration  vs  Time’/ 
end;  ’ •  d: sp 1 ay_p 1 o t s  •> 


procedure  integrate  (method  :  char;  step,  ma>.t;me  :  reaii  var  i  :  integer); 


<  *  ♦ ) 

(*  Abstract:  This  subroutine  carries  out  the  numerical  integration  •> 

<•  process.  If  method  •  then  the  linear  acceleration  method  is  •) 

(*  used.  Otherwise,  the  constant-average-eccelerat ion  method  is  used.  *> 

(•  step  specifies  the  time  step  to  be  used  and  maxtime  specifies  the  •) 

<«  length  of  time  over  which  to  integrate,  t  returns  the  total  number  •) 

(*  of  points  that  were  integrated  over.  In  addition  to  using  the  •  > 

<*  matrix  routines  provided  by  mopbox  ,  integrate  calls  various  user-  •> 

<*  defined  routines  to  get  information  about  the  particular  problem  *) 

(*  at  hand.  The  following  routines  are  called  from  the  file  •) 

<•  ’USER. TEXT’:  •> 

<*  get_ndofs  -  returns  the  number  of  degrees  of  freedom  *) 

<*  get_inl t_values  -  returns  Initial  displs,  velocs,  and  accels  •) 

<*  get_mass  -  returns  the  mass  matrix  •> 

(•  get_st i f f ness  -  returns  the  etiffness  matrix  ») 

<*  get_damping  -  returns  the  damping  matrix  •) 

<*  get_force  -  returns  the  force  vector  *) 

(•  w) 

. . * . * . * . . 


type 

int_const  -  arrayt0..7)  of  reali 
var 

T  :  real) 

Dp,Up,Ap  :  vector) 

Dc ,Uc  .Ac  :  vector i 

mass  :  matrix) 
damping  :  matrix) 
stiffness,  Keff  :  matrix) 

F,  Faff  .*  vector) 
a  :  int_const) 

vl .v2.v3.v4.v5  !  vector) 
ml  ,m2  :  matrix ) 

count  !  integer) 
npts,  nskip  :  integer) 

procedure  save_data<T  s  real)  0.  U,  A  :  vector i  i  :  integer)) 

(*  saves  the  current  data  'in  the  global  arrays  *) 
var 

dof  :  integeri 
begin 

t imal i  1  :■  T  *  1000.0) 
for  dof  :*•  1  to  ndofs  do 
begin 

datat dof  1. displt 1 1  :•  Dtdoflt 
datat do f 1 . veloct i  1  :■  VCdofli 
datat dof 1 .accell i  ]  :■  Atdof] 
end 

end)  (•  save_data  *> 

procedure  get_conatants  (method  :  char>  step  :  reali  var  a  :  int_const)i 

(*  returns  the  integration  constants  for  newmark  integration  •) 

var 

alpha  ,  delta  :  real  i 
begin 

if  (method  -  'I')  then  <*  linear  acce lera  t ion  method  *  ) 
alpha  : ■  1.0  /  6.0 

else  (•  constant  jvorage  acceleration  »> 
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alpha  :«  0.25i 
delta  :■  0 . 5 i 

a(0]  1.0  l  (alpha  •  3or(step>>» 

alii  oelta  /  (alpha  *  step)i 

a!2 1  :«  1.0/  (alpha  •  step)) 

.13)  j-  1 .0  /  (2.0  *  alpha)  -  1 .0t 
s[  4 )  :■  delta  /  alpha  -  1.0i 
aCS)  !•  step  /  2.0  •  (delta  /  alpha  -  2.0)i 
at 6)  :«  step  *  (1.0  -  delta)! 
aC7)  delta  •  step 

end»  (•  get_constants  «> 

begin  <•  integrate  •) 


ndofs  :■  get_ndofsi  <* 

if  set_si2e(  ndofs  )  them  (• 

get_ini t_values< Cp ,  Op,  Ap)» 
get_constants( method .  step,  a)i 
count  :*  1 j 

i  !■  01 

T  i m  0.0i 


get  the  numb  r  of  degrees  of  freedom  •> 
set  size  of  matrix  operations  in  mopbox  •  ) 
(»  get  the  initial  values  •> 

(•  get  newmark  Integration  constants  •> 
(»  counts  each  time  step  taken  •) 

(«  the  current  number  of  plotting  pts  *  ) 
(*  start  counting  at  time  T  -  0  •) 


npts  :»  trunct maxt ime  /  step)  +  I j 
nskip  :»  trunc(npts  /  MAXPTS >  +  li 


while  <T  <■  maxtime)  do 
begin 

if  (< count-1)  mod  nskip)  •  0  then 
begin 

i  J  *  i  +  1  « 

save_data(T,  Op.  Op,  Ap ,  i)  (*  for  output  and  plotting  •  ) 
endt 


(•«••  g8t  the  mass,  damping,  stiffness,  and  force  matrices  •••* > 

get_mass(Dp,  Op,  mass/i 
get_dampmg(  Dp  ,  Op,  damping)! 
get_s t i f f ness( Dp ,  Op,  stiffness)) 
get_force( Dp ,  Op,  T,  F>i 

<♦*•*  form  effective  stiffness  matrix  ••««) 

<*  Keff  :*  stiffness  +  a[0)»mass  +  a[))*damping  •) 

scalem( mass ,a(0] ,ml  )i 
scaiem(  damp  mg  ,at  1  3  ,m2  )  i 
addr.f  ml  ,m3  ,ml  )  i 
addm( s 1 1 f f ness  ,ml  ,Keff)( 

( form  effective  load  vector  *«**) 

(*  Ferf  :■  F  +  mass*( at  0  ]*Dp  +  atZl’Op  +  aC3)»Ap> 

+  damping*( aC 1  )«Dp  +  a(4]*0p  +  aC53»Ap>  •) 

scalev(Dp,a[03  , v J  >i 
sea  le  v<  Op  ,aC  2 ) ,  v2  )  j 
scalev!  Ap  ,a[  .3  )  ,  v3  ) ! 
addv ( v 1  ,  vZ  ,  v l  )  f 

addv( v 1  ,  v3  ,  v 1  )  i 
mul t ( ma ss  ,  v 1  ,  v4  ) » 


54 


scalevt Op ,a( 1 1 , vl > t 
scalevt Up ,at 43 ,v? >i 
scalevt Ap ,a(53 ,v3  >i 
addv<  v I , v2 , v  I  ) i 
addvt vl  ,v3,vt  >i 
multtdamping.vf  tvB>i 

addv(F,v4,Feff  >i 
addvtFaf f ,v5 ,Fef f )» 

(••••  aolva  for  displacements  •*••) 

<*  Oc  r-  Kef f*-l  *  Faff  *> 

if  net  sol vet  Kef f ,Fof f ,0c >  then 

write:./ ’can’ ' t  solve  equation  s/stem  (zero  in  diagonal >'>t 

<«*•*  determine  accelerations  •#•«) 

(*  Ac  ! “  al03*<Oc  -  Dp)  ~  at 23«Up  -  at33*Ap  *) 

subv(Dc  ,Dp ,vl  >* 
scalav( vl ,a[03 ,vl )t 
sea  lev(  Up  ,at  21  ,v2  >♦ 
scalev( Ap ,al 31 ,v3  )i 
subv( vl  ,  v2  ,  V I  )  | 
subv<  vl  ,v3.Ac  )i 

(««*•  determine  velocities  «*••) 

(*  Uc  S“  Up  +  e(63*Ap  ♦  a£7]«Ac  *> 

scalevt  Ap ,e(Gl  ,vl  >» 
scalevt Ac ,a£7J ,vZ  )» 

addv<  Up  ,  vl  ,  vl  )  i 
addv<  vl ,v2.Uc)i 

<•*«*  reassign  current  values  to  previous  variables  •*#• ) 

copyvt Dc  ,Dp > j  (•  Dc  ->  Dp  •> 

copy v(  Dc ,Up ) i 
copyvt Ac ,Ap ) i 

<*•««  increment  time  by  one  time  step  ***♦ ) 

count  >  *  count  +  I  » 

T  : «  t  count  -  t  )  *  step i 
end 

endi  ( »  integrate  •  ) 


begin  ( *  newmark  • > 

wri taint '  NEWMARK  'It 

writelnt  '  Numerical  Integration  Program  ’)i 

repast 

wri tel n i 

wn t el n< ' Solut i on  method 

writelnt'  I)  Linear  Acceleration  Method’)! 
writelnt'  2)  Cone tant -Averaga-Acceierat ion  Method'  )i 
writet  'Enter  the  number  desired  !  ’>! 
readt  ch  )  t 
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writalni 

until  <ch  in  [T,'2’l)| 

writei 'Enter  ine  step  size  (mSec):  *  > * 

readlnt a tep_size ) t 

fltap^size  :*  0.001  •  stop_aizei 

wri  te-t  ’  Enter  the  length  of  tine  (msec):  ’)» 

readln<  next ime  ) i 

maxtime  I “  0.001  •  maxtinei 

lntegratet ch ,  atep_»ize,  maxtime,  npta)i 

repeat 

wrlte< 'Output  data  to  uhat  file  ( <enter>  ignores)  :  ' > i 
readln< ? i lenane > i 

until  ( filename  ■  ' '  >  or  open_out fi le< fi lenane > i 

if  < filename  <>  '')  then 
begin 

writelnt 'Writing  to  file  ',  filename)* 
outpu t_data( npt a  > i 
closet  out  file ,  'save') 
end  i 

wri taint 'Enter  dof  number  to  view  t*'q*'  to  quit)')* 

repeat 

wr itet ' >  '  )| 
raadt  ch  >i 
writalni 

if  (ctoi(ch)  in  Cl..ndofa))  then 
diapla/.plotat ctoit ch  ) ,  npts)i 
until  (ch  In  C 'q'  ,  ’ 0 ’ 3  > I 

end.  <•  newnerK  ♦> 
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module  mopboxi 

. . . . . . .  ) 

<  •  *  ) 

(*  flflFBOX  •  > 

<•  Matrix  Operations  toolBOX  •) 

<•  •) 

<•  Author:  John  Baugh,  Battelle-Northwest  Laboratories  •> 

(•  Data:  November  1,  1985  *) 

( •  *  ) 

,<•  Abstract:  This  module  provides  low-level  operations  to  manipulate  *> 

_<*  matrix  and  vector  entitles.  In  order  to  use  the  functions,  the  *5 

<•  calling  program  must  first  call  set_size  to  define  the  order  of  •) 

<•  the  matrices.  For  instance,  the  call  *) 

(  *  if  set_eize(8)  then  .  ..*  ») 

<•  from  the  main  program  results  In  all  matrices  being  defined  as  •> 

(•  8x8  and  all  vectors  as  8x1.  It  should  be  noted  that  attempts  •) 

<•  to  set  the  size  larger  than  the  maximum  declared  size  (10)  •> 

<*  results  in  set_eize  returning  false  (and  the  size  remaining  O 

<*  unchanged).  *> 

<  •  ♦  ) 

<•  The  calling  program  must  Import  this  module  before  it  can  access  *) 

<•  any  of  its  functions,  which  may  be  done  using  the  ‘import*  compiler  •) 

( *  directive.  The  nodule  object  code  must  be  online  during  compilation  *) 

<•  of  the  calling  program  (use  the  'search*  compiler  directive  or  *) 

<•  include  this  module  in  the  system  library  —  see  chapter  I  of  the  HP  •) 

<•  Pascal  3.0  Procedure  Library  for  more  details).  Before  executing  •) 

<*  your  program,  permanently  load  this  nodule  into  memory  using  the  *) 

<•  p-load  command  (see  page  120  of  the  HP  Pascal  3.0  User's  Guide  for  •) 

<•  more  information).  •> 

<  •  •  ) 

( •••••••••••«•• •*••«••••••••••••••»••••••••«••••• •«•«•••••*  ) 

export 

const 

MAX  »  1 0 1  <•  maximum  size  of  matrices  and  vectors  •> 

type 

vector  ■  array  U..MAX1  of  real* 
matrix  ■  array  [I., MAX,  I..MAXI  of  reali 

function  set_slze  (size  i  Integer)  i  boolean! 

(•  set  the  order  for  ALL  matr ix/vector  operations  in  the  mopbox  module  •) 

function  solve  (ml  :  matrix)  vl  :  vector)  var  v2  >  vector)  :  boolean* 

(•  solve  up  to  10  simultaneous  linear  algebraic  equations  *) 

procedure  mult  (ml  s  matrlxi  vl  :  vector)  var  v2  :  vector)) 

(«  multiply  the  matrix  ml  by  vector  vis  v2*mt«v1  •> 

procedure  zerom  (var  ml  :  matrix)) 

(•  zero  each  element  in  matrix  ml  •) 

procedure  2erov  (var  vl  :  vector)) 

(*  zero  each  element  in  vector  vl  •) 

procedure  copym  (ml  :  matrix)  var  m2  :  matrix)) 

(•  copy  matrix  ml  into  matrix  m2  *) 

procedure  copyv  (vl  :  vector)  var  v2  :  vector)* 

(•  copy  vector  vl  into  vector  v2  •) 
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procedure  acalem  (ml  :  Matrix i  fl  :  reali  var  m2  t  Matrix )i 
'ft  •  scale  each  element  in  matrix  ml  by  fl,  puttinQ  it  into  matrix  m2  *> 

procedure  scalev  (vl  :  vactori  fl  :  real*  var  v2  t  vector >t 

<•  scale  each  element  in  vector  vl  by  fl ,  puttinQ  it  into  vector  v2  *> 

-procedure  addm  <ml ,  m2  :  matrix i  var  m3  i  matrix)! 

(*  add  matrix  m2  to  matrix  ml:  m3«mt+m2  *> 

procedure  addv  <vl ,  v2  :  vector i  var  v3  :  vector >i 
<*  add  vector  v2  to  vector  vis  v3«v1+v2  *) 

.procedure  subm  (ml,  m2  :  matrix!  var  m3  :  matrix)! 

(•  subtraot  matrix  m2  from  matrix  ml:  m3“ml-m2  ♦) 

procedure  subv  (vl  ,  v2  :  vector!  var  v3  :  vector)! 

(•  subtract  vector  v2  from  vector  vl :  v3*v1-v2  »> 

implement 

var 

n  :  lnteoari  (•  order  of  matrix/vector  operations  •> 


function  set_size  (size  :  mteoer)  :  boolean! 

(•  set  the  order  for  ALL  matr ix/vector  operations  in  the  mopbox  module  *) 
beQin 

if  (size  <■  MAX )  and  (size  >  0)  then 
begin 

set_slze  :■  truei 
n  :•  slzei 
end 

else 

set_aize  :■  falsei 
end i  <  *  set_size  • ) 


function  solve  (ml  :  matrlxi  vl  :  vectori  var  vZ  ?  vector)  :  boolean! 

<••*♦* . . * . .  ) 

<  •  •  > 

<*  Abstract:  This  routine  solves  up  to  10  simultaneous  linear  elQebralc  •) 
(•  equations  by  Gauss-Jordan  Elimination.  If  all  of  the  diagonal  *> 

(•  terms  in  the  coefficient  matrix  are  non-zero,  solve  returns  true.  *> 

<♦  In  general,  for  an  accurate  solution,  the  diagonal  terms  should  *) 

<*  dominate  the  coefficient  matrix,  ml.  *) 

( »  •  ) 

(*  input  parameters  -  #> 

(•  ml  :  coefficient  matrix  •) 

(•  vl  :  right-hand  vector  •) 

(*  output  parameters  -  •> 

(•  v2  :  solution  vector  •) 

<•  solve  :  returns  trus  if  all  terms  in  the  diagonal  *> 

(•  of  the  coefficient  matrix  are  non-zero  *> 

( •  •  ) 

. . . . - . * . . 

label 

I  ) 

var 

1  ,  J  ,  k  :  integer! 


hold  :  reali 
begin 

solve  :■  true* 

for  i  :•  1  to  n  do  (*  loop  over  each  of  the  columns  of  ml  •) 

begin 

for  k  1  to  n  do  <•  loop  over  each  row  except  pivot  row  *) 
begin 

if  ( k  <>  i  )  then 
begin 

if  (nl[l,ll  <>  0.0)  then  <•  be  sure  diagonal  is  non-zero  •) 
begin 

hold  «•  -n1tk.il/m1ti, ill 

for  j  :•  I  to  n  do  (•  loop  over  each  element  In  a  row  •) 

begin 

mltk.jl  :■  mltk.jl  +  hold  •  mlti.jl! 
if  <J  ■  i>  then 
mltk.jl  :■  0.0 

endi 

vltk]  :»  vlCkl  +  hold  •  vltil 
end 

else  <•  found  a  zero  diagonal  tern  In  coefficient  matrix,  ml  •> 
begin 

solve  :■  falsai 
goto  1 i 
end 
end 
and  i 

hold  i-  ml t i  ,1 ] » 

for  J  i-  1  to  n  do  <*  this  loop  is  for  each  element  in  •> 

mlti.jl  :•  mlti.jl  /  hold!  <•  the  pivot  row  •) 

mltt.il  :•  1.0i 
vltil  vltil  /  hold 
endi 

for  l  :•  I  to  n  do  <•  copy  solution  into  solution  into  •) 

v2tll  :*  vltili  (•  the  solution  vector  •) 

1  : 

endi  < •  solve  • ) 

procedure  mult  (ml  i  matrix!  vl  5  vector!  ver  v2  »  vector)! 

<•  multiply  the  matrix  ml  by  vector  vl:  v2anl*v1  •) 
ver 

i  ,  j  :  integer! 
begin 

for  l  :«  I  to  n  do 
begin 

v2til  0.0i 
for  j  :■  I  to  n  do 
v2til  :•  v2C i 3  +  mlti.jl  •  vltjl 

end 

endi  (»  mult  •) 


procedure  zerom  (ver  ml  :  matrix)! 

<•  zero  each  element  in  matrix  ml  •) 
ver 

1  ,  j  :  integer! 
begin 

for  i  : ■  I  to  n  do 
for  j  : ■  I  to  n  do 
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ml  I i  ,  j  1  :•  0.0 
endi  <«  zeron  •) 


procedure  zerov  (var  vl  :  vector )i 
<•  zero  eech  element  in  vector  vt  •) 
ver 

i  :  integer! 
begin 

for  i  !■  1  to  n  do 
v  i  [  i  ]  :<*  0.0 
end i  (•  zerov  *) 


procedure  copym  (ml  i  matrix t  var  m2  s  matrix >i 

(*  copy  matrix  ml  into  matrix  m2  •> 

var 

1  ,  J  :  integer i 
begin 

for  i  :*  I  to  n  do 
for  J  :■  I  to  n  do 
m2C i  ,  j  1  :•  ml  [  i  ,  J  ] 
endi  <•  copym  «) 


procedure  copyv  (vl  :  vector  i  var  v2  :  vector  )| 

( *  copy  vector  vl  into  vector  v2  ♦  ) 

var 

1  :  integer) 
beQin 

for  i  :»  I  to  n  do 
v2Cil  :•  v I C i J 
endi  ( •  copyv  *  ) 


procedure  scalar  (ml  :  matrix)  fl  :  real!  var  m2  :  matrix)! 

(♦  scale  each  element  in  matrix  ml  by  fl .  putting  it  into  matrix  m2  •) 
var 

1 .  j  :  integer! 

br  ->in 

for  i  :*  1  to  n  do 
for  j  : ■  I  to  n  do 

m2! l  ,j  1  :■  f 1  •  ml ( i  ,J  ] 
end i  < •  scalar  •  ) 


procedure  scalev  (vl  :  vector!  fl  :  real!  var  v2  :  vector)! 

(•  scale  each  element  in  vector  vl  by  fl,  putting  it  into  vector  v2  •) 
var 

i  :  integer! 
begin 

for  l  :**  1  to  n  do 
v2C ll  :■  fl  •  vlt i J 
endi  < •  scalev  •  ) 


procedure  addm  (ml ,  m2  :  matrix!  var  m3  :  matrix)! 

<»  add  matrix  m2  to  matrix  ml:  m3*m!+m2  *) 

var 

l  ,  J  :  integer) 
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for  i  i»  1  to  n  do 
for  j  t-  1  to  n  do 
m3ti,j3  !•  mlti.j]  ♦  m2(  l  ,J 3 

endi  <*  addn  •) 


procedure  addv  <  v l ,  v2  :  vector i  ver  v3  :  vector )» 
(•  add  vector  v2  to  vector  v 1 t  v3*vl+v2  *> 

J  var 

i  »  integer i 
begin 

for  i  :•  t  to  n  do 

v3C 1 J  :•  vltil  +  v2ti] 

end i  <  *  addv  «  > 


procedure  subm  (ml,  n2  :  matrix)  var  n3  t  matrix >i 
<«  subtract  matrix  m2  from  matrix  ml «  m3«m1-m2  *) 
var 

1,  J  :  integer i 
begin 

for  i  »■  1  to  n  do 
for  J  :•  I  to  n  do 
m3ti,J]  :•  mICljl  -  m2ti,J] 
end*  <•  subm  •) 


procedure  subv  (vl,  v2  :  vectori  var  v3  :  vector)| 
(•  subtract  vector  v2  from  vector  vli  v3»v1-v2  *) 
var 

1  J  integer i 
begin 

for  l  :»  1  to  n  do 

v3C 1 1  :•  vl  C 1 3  -  v2t  1 3 
endi  <•  subv  *) 

end.  <•  mopbox  •) 
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nodule  plotpeci 
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PLOTP/IC 
Graphic  Plotting  Package 

Author:  John  Baugh,  Bat  tel le-Northwest  Laboratories 

Date:  August  12,  1 985 

Abstract:  This  is  a  library  for  plotting  data  in  the  form  of  line 
graphs.  The  only  information  needed  by  the  plotting  routine  is 
the  arrays  of  points,  the  number  of  points,  an  index  for  the  line 
style  (e.g.  solid*!  ,  dotted-7,  etc.),  and  a  title  which,  when  null, 
causes  the  current  plot  to  be  displayed  over  the  previous  one. 

Plots  are  automatically  scaled  for  convenience.  A  typical  example 
is  illustrated  below: 

program  example  < input,  output): 
esearch  ' VOLUME : PLOTPAC ' • 
import  plotpaei 
const 

SOLID  -  1 
DASHED  -  2 
var 

x,  y,  z  :  plot_arrayi 
keyboard  :  texti 
ch  :  chari 
begin 

resat( keyboard  ,  'console:* >i 

<  define  x,  y,  and  z  arrays 

from  1  to  numpts  here  > 


and  z  versus 


t’  >i 


if  inltgraphics  then 

begin 

plottx,  y,  numpts,  SOLID,  ’y 
plot(x,  z,  numpts,  DASHED,  * 
read( keyboard ,  ch>» 
termgraphics 
end 
end. 


As  in  the  example,  the  calling  program  must  initialize  the  graphics 
display  by  calling  "inltgraphics"  before  plotting  data  with 
procedure  "plot".  After  plotting  the  data,  a  “read"  or  "readln" 
statement  located  before  a  call  to  “ termgraphics *  keeps  the  plot 
displayed  until  the  user  gives  an  appropriate  response.  In  this 
case,  typing  any  key  will  remove  the  plot  and  return  the  display  to 
alpha. 

The  calling  program  must  import  this  module  before  it  can  access 
any  of  its  functions,  which  may  be  done  using  the  "import"  compiler 
directive.  The  module  object  code  must  be  online  during  compilation 
of  the  calling  program  (use  the  “search"  compiler  directive  or 
include  this  module  in  the  system  library  —  see  chapter  1  of  the  HP 
Pascal  3.0  Procedure  Library  for  more  details!.  Before  executing 
your  program,  permanently  loed  this  module  into  memory  using  the 
p-load  command  (see  page  120  of  the  HP  Pascal  3.0  User’s  Guide  for 
more  information). 
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import  dgl_llb*  <*  access  graphics  library  *) 

export 

const 

WAXPTS  -  150* 
type 

plot_string  ■  stringC803» 

plot_array  •  array  C1..MAXPTS3  of  real * 

function  initgraphics  :  boolean! 

(*  initialize  tha  device-independent  graphics  library  and  display  *) 
procedure  termgraphics i 

<•  terminate  the  device-independent  graphics  library  and  display  •) 

procedure  plot  ( var  x,  y  :  plot_arrayi  n  J  integer*  linetype  ’•  integer* 
title  :  plot_string >i 
<•  plot  the  input  arrays  *) 

implement 

const 

ASPECT  -  3.0*  <•  aspect  ratio  for  all  characters  •) 

type 

modetype  -  (alpha,  graphics)* 
roundtype  ■  (up,  down,  naar)i 
var 

initstate  :  boolean*  <•  true  if  graphics  initialized  •) 

xmin,  xmax  s  real*  <*  min  and  max  values  in  x  array  «) 

ymin,  ymax  :  real*  <•  min  and  max  values  in  y  array  •) 

. . . . . .  math  ROUTINES  . . . . 

function  log  <x  ;  real)  s  reali 
<•  return  log  base  10  of  x  *) 
begin 

log  ln( x  )  /  ln< 10.0  ) 
end*  <  *  log  •  ) 


function  power  (base,  exponent  :  real)  '•  real* 
(•  return  the  value  of  base  *  exponent  •  ) 
begin 

if  (base  ■  0.0)  then 
power  :■  0.0 
else 

1  f  ( base  >  0.0)  t hen 

power  :•  exp(ln(base)  •  exponent) 
else 

if  (base  <  0.0)  then 

power  :•  -  exp( ln< -base >  •  exponent) 
end*  ( *  power  »  ) 


procedure  scinotation  (x  :  real*  var  mantissa,  exponent  :  real)) 
(»  return  x  in  the  form  of  scientific  notation  •> 
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bag  in 

if  (x  •  0.0)  then 
begin 

mantissa  :■  0.0) 
exponent  :■  0.0 
end 

^elee 

if  (x  <  0.0)  then 
begin 

mantissa  «■  -power<10,0,  log(-x)  -  trunct log( -x > > ) i 
exponent  :•  trunct log< -x > > 
end 
else 

if  (x  >  0.0)  then 
begin 

mantissa  :«  power(10.0,  log<x>  -  trunc( log( x  )  > )  i 
exponent  :•  trunc< log< x ) ) 
endt 

if  < abs< mant issa  )  <  1.0)  then  (•  normalize  the  result  •) 

begin 

mantissa  :■  mantissa  *  10. 0t 
exponent  : ■  exponent  -  1.0 
end 

end)  <*  scinotation  *) 


function  absmax  <x,  y  :  real)  :  real  I 

<•  return  the  maximum  of  the  absolute  values  of  x  and  y  •) 
begin 

if  (abs(x)  >  abs(y))  then 
absmax  : *  abs<  x  ) 
else 

absmax  :■  abs< y  ) 
end i  <  *  absmax  *  ) 


function  roundZ  (n,  m  :  real)  mode  :  roundtype >  •  real) 

<*  round  n  to  the  nearest  m,  according  to  mode  •> 
const 

ROUNDOFF  ■  IE-100)  (*  roundoff  error  fudge  factor  •  > 


var 


rounded  :  reali 
negative  :  boolean) 
begin 

negative  :■  (n  <  0.0  )| 
if  negative  then 
begin 

n  :■  abs( n  )  i 
if  mode  ■  up  then 
mode  :•  down 
else 

if  mode  •  down  then 
mode  :■  up 


<*  temporary  holding  area  •) 
<*  flag:  “is  it  negative?"  •) 

<»  is  the  number  negative?  •) 


(•  work  with  a  positive  number 


end) 

case  mode  of 

down  :  rounded  :*•  truncfn  /  m)  *  mi 
up  :  beg.n 

rounded  :»  n  /  m  i 

if  abs( rounded  -  round< rounded  )  )  >  ROUNDOFF  then 
rounded  :•  ( trunc< rounded )  +  1.0)  •  m 


•  ) 
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•  Isa 

rounded  :■  trunc< rounded )  *  n 
end  i 

near  :  rounded  :•  trunc<n  /  n  ♦  m  •  0.5)  •  m 
end  i 

if  negative  then  (*  reinstate  the  sign  *) 

rounded  :•  -roundedi 
round2  !■  rounded 
end »  <•  round2  •> 

. . . . .  MACHINE-OEPENOENT  GRAPHICS  ROUTINES  •••**•*•* 

function  setnode  (node  :  nodetype)  :  boolean) 

<»  set  the  display  mode  to  alpha  or  graphics  •> 
const 

AlPHAOISP  ■  1 05 1 1  (*  alpha  display  address  •> 

6RAPHDISP  -  10501  <*  graphics  display  address  •  > 

var 

error  :  integer i  (•  output_esc  error  (0  if  ok)  •) 

on,  off  s  integer i  <*  on/off  switches  «) 

x  :  real;  <•  ignored  •) 

begin 

setnode  )■  true) 
on  :•  It 
off  :■  0i 

if  (node  •  alpha)  then 
begin 

output_esc< GRAPHOISP ,  1,0,  off,  x,  error))  (•  graphics  off  •) 
if  (error  <>  0)  then 
setnode  )•  falsei 

output_esc( ALPHAOISP ,  1,  0,  on,  x,  error  )|  (•  alpha  on  *) 

if  (error  O  0)  then 
setnode  ;•  false 

end 

else 

begin 

output_esc( ALPHAOISP ,  1,0,  off,  x,  error))  (•  alpha  off  *> 
if  (error  <>  0)  then 
setnode  :■  false) 

output_esc( GRAPHDISP ,  1,  0,  on,  x,  error))  (•  graphics  on  •) 
if  (error  <>  0)  then 
setnode  :■  false 

end 

end)  <  *  setnode  * ) 


function  initgraphics  :  boolean) 

<•  initialize  the  device-independent  graphics  library  and. display  •> 
const 

XPIXEL  ■  400i  (•  pixels  in  x  dir  <  HP  98 1 6  >  •) 

YPIXEL  -  300)  <*  pixels  in  y  dir  (HP9816)  *> 

OEUICE  -  3)  <«  address  of  screen  *) 

CONTROL  "0i  <•  ignored  by  screen  *) 

var 

error  :  integeri  (•  display  imt  error  (0  if  ok)  •) 

begin 

graphics_ini t i 

dlsplay_init(OEUICE ,  CONTROL,  error)) 
if  (error  ■  0)  and  setnode! graphics )  then 
begin 


I 


lnitstate  :«  true? 

•et_aspect< XPIXEL  -  I  ,  YPIXEL  -  I  )  !•  use  entire  screen  ♦  > 

end 

else 

lnitstate  :■  falsei 
lnitgraphics  :»  lnitstate 
endt  <•  lnitgraphics  *> 


procedure  tarmgraphicsi 

<•  terminate  the  device-independent  graphics  library  and  display  •) 
begin 

if  lnitstate  then 
begin 

If  setmode<alpha  )  then 
begin 

graphics_termi  <•  terminate  graphics  device  *) 

lnitstate  :•  false  <*  uninitialize  graphics  *) 

end 
end 

endi  <*  termgraphlcs  *> 

<••••••••••••«•••«••••«•«••*«  PLOTTING  functions 

procedure  setcharsize  (width  :  real)) 

(*  set  the  character  size  using  ASPECT  as  the  aspect  ratio  •) 
begin 

set_char_size( width ,  ASPECT  *  width) 
endi  <*  setcharsize  ») 


procedure  get_format  (x  :  reaii  var  numlen,  fraclen  :  integer  )i 
(*  determine  the  required  string  length  (total  and  fraction)  of  x  •  > 
begin 

i f  abs<  x  )  >  0  then 

numlen  :•  trunc( log! abs( x  ) )  +  I) 
else 

numlen  :■  0t 
fraclen  :»  4  -  numleni 
if  numlen  <  I  then 
numlen  : •  I i 
if  fraclen  <■  0  then 
begin 

fraclen  '• "  0t 

numlen  :*  numlen  +  fraclen  +  li 
end 
else 

numlen  :»  numlen  ♦  fraclen  2» 
endi  (*  get_format  •) 


procedure  writenum  (x  :  reali  numlen,  fraclen  :  integer)) 

<  *  write  x  on  the  graphics  display  at  the  current  location  •) 

var 

slen  :  integeri 
s  :  plot_strmgi 
begin 

strwrite(s,  I,  slen,  *: numlen: free len >; 
setstrlen<3,  slen  -  ))i 
gt ex t ( s  ) 
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ana i  < *  wr 1 tenum  *  ) 


ppocadura  writeatr  (a  t  plot_atnng ) i 

(♦  write  etring  a,  centered  and  at  the  top  of  the  acraen  •) 
conat 

WIDTH  ■  0.04t 
begin 

aetcharaize<WIOTH )i 

movet "8trlen( a >  •  WIDTH  /  2.0,  0.9)i 

gtext<  a ) 

endi  ( *  writeatr  • ) 


procedure  draubox  (x,  y  :  real)! 

(•  draw  a  box  through  pointa  <x,y)  and  <-x,-y>  •) 
begin 

move<  x  ,  y  )  i 
line(  x  , -y  > i 
1 ine( ~x  ,-y  )  i 
1  ine< "x  ,  y  )  i 
1 ine<  x ,  y  ) 
end i  ( •  drawbox  •  ) 


procedure  findmaxmin  (var  a  :  plot_arrayi  n  :  integeri  var  max,  min  :  real)! 

<*  determine  the  max  and  min  valued  of  array  "a"  relative  to  zero  •) 

var 

i  :  integeri 
begin 

max  :  ■  0 . 0 1 

min  :■  0.0t 

for  i  I  to  n  do 

begin 

if  at i J  >*  max  then 
max  : ■  aC i 1 t 
if  a[ i  1  <■  min  then 
min  : *  at  1 1 1 

end! 

if  (max  *  0.0)  and  (min  ■  0.0)  then 
max  : ■  1  . 0 1 
endi  (*  findmaxmln  •) 


procedure  f indapacing< var  max,  min,  spacing  :  realh 

(•  determine  the  apacing  and  adjuat  the  max  and  min  vaiues  appropriately  *> 
var 

mantissa,  exponent  :  reali  (*  mantlaaa  and  exponent  of  abamax(x)  •) 
a  :  integeri  (*  a  &  b  are  uaed  to  determine  apacing  •) 

b  :  reali 
begin 

sc i notat i on< absmax ( max ,  min),  mantlaaa,  exponent  )| 
a  :•  trunc( mant l asa  )  +  I l 
b  :*  power(10.0,  exponent)! 
case  a  of 


1 

spacing 

J  a 

0.20 

# 

b  i 

2 

apacing 

•  m 

0 . 25 

* 

b  i 

3,4 

spacing 

;  m 

0.50 

• 

b  i 

5  ,6  ,7,8 

apacing 

;  « 

1  .00 

• 

b  i 

9,10 

spacing 

*  m 

2.00 

• 

b 

e-.di 

max  round2<nax,  spacing,  up)( 
run  :«  rounfl‘.((nn,  spacing,  down) 
andi  '•  flndipaclng  ») 


function  glx  (xglobel  s  real)  :  realt 

<*  convert  global  x  coordinates  to  local  x  coordinates  *) 
begin 

glx  2.0  *  (xglobal  -  xrnin)  /  < xnax  -  xrnin)  -  1.0 
endi  (*  glx  ♦) 


function  gly  (yglobal  :  real)  :  reali 

<*  convert  global  y  coordinates  to  local  y  coordinates  «) 
begin 

gly  :■  2.0  *  (yglobal  -  ymin)  /  (ymax  -  ynin)  -  )  .0 
endi  <•  gly  •) 


procedure  dra  xes i 

(*  draw  the  x  and  y  axes  *) 

begin 

novel  glx<  xrun  ) ,  p'.(0.0))i 
line!  glx<  xnax  > ,  gly<0.0))» 
novel glx < 0 . 0 >  gly< ynin  )  )| 
line<  glx<  0.0  :  ,  gly<ymax>> 
er.di  ( *  drawaxes  •  ) 


function  i n_range<  x ,  y  :  real)  :  booleani 

1 •  determine  whether  or  not  point  (x,y)  lies  inside  the  display  area  * 
begin 

if  (x  <>  xnax)  and  (x  >■  xnin)  and  (y  <■  ymax)  and  <y  ynin)  then 
in_range  :■  true 
else 

in_range  ;■  false 
end i  < *  in_ranoe  •  ) 


procedure  drawlines  ( var  x ,  y  :  plot.^arrayi  n  :  integer)) 

<•  draw  lines  connecting  the  points  in  the  x  and  y  arrays  •) 
var 

l  :  lntegeri 
count  :  integer) 
begin 

count  < “  I ; 

while  (not  in_range< x[ count  ] ,  ylcount]))  and  (count  <  n)  do 
count  J "  count  +  1  i 
novel glx( xtcount 1 ) ,  gly'yCcount]  )  )i 
1  f  count  n  t hen 

for  i  : *  (count  +  1  )  to  n  do 
if  in_range< x( l 1 ,  ytil)  then 
l.ne(gix( xC i ]  ) ,  gly(ytil)) 
end )  ( »  drawl ines  * ) 


procedure  labelxaxis  (<spacmg  :  real)) 

<♦  put  tick  narks  and  scale  numbers  on  the  x  axis  *> 
const 
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WIDTH  -  0.03i  («  width  of  a  character  (loeal)  *) 

SEMI  •  0.02)  <•  half  length  of  a  major  tick  nark  *) 


var 

xg  :  real i  (  * 

xl ,  y 1  s  real t  <  • 

nunlen,  fraclen  :  integer)  <« 
height  :  reali  <* 

begin 

setcharsize( WIOTH) t 
height  :■  ASPECT  •  WIDTH) 


get_fornat(absnax<xnax  ,  xnin), 
yl  :■  gly<  0.0  )i 
xg  :■  xeuni 


global  x  coordinate  •) 
local  x  and  y  coordinates  *) 
number  and  fraction  lengths  •) 
character  height  *) 


nunlen,  fraclen)i 


while  <xg  <»  (xmax  +  0.S  *  xspacing))  do 
begin 

if  abs(xg)  <  (0.S  •  xspacing)  then 

if  (ymin  <  0.0)  and  (ynax  >  0.0)  then 
xg  :■  xg  +  xspacingi 
xl  glx( xg  > i 
move<  xl  ,  yl  -  SEMI  )| 


< •  x  is  zero  *  ) 

(•  y  axis  crosses  x  axis  •) 
(•  skip  to  the  next  label  •> 
(♦  convert  x  from  global  to  local  •) 
(•  move  to  the  proper  location  «) 


1  ine( x 1  ,  y 1  ♦  SEMI  ) i 


<•  and  draw  the  tick  nark  •  ) 


i f  <  ynax  O  0 . G > 
move! xl  -  nu' 

else 

novel xl  -  nunlen  * 


<«  x  axis  at  top  of  the  screen  •) 
WIDTH  /  2.0,  yl  +  0.S  •  height) 

(•  x  axis  is  somewhere  below  top  •) 
WIDTH  /  2.0,  yl  -  1.5  •  height )i 


writenuml  xg ,  nunlen,  fraclen)i 
xg  :•  xg  +•  xspacing 
end 

endi  I •  labelxexis  *  ) 


procedure  labelyaxis  (yspaclng  :  real); 

I*  put  tick  narks  and  scale  numbers  on  the  y  axis  *) 
const 

WIOTH  ■  0.03)  <*  width  of  a  character  (local)  •> 

SEMI  -  0.02i  <*  half  length  of  a  major  tick  nark  •) 

var 


yo 

:  rea 

1) 

(  * 

global  y  coordinate 

•  ) 

xl  , 

yl  : 

real  i 

<  * 

local  x  and  y  coordinates  •) 

nunlen , 

fraclen  :  integer i 

<  • 

number  and  fraction 

lengths  • ) 

height  : 

real  i 

(• 

character  height  •) 

begin 

set  char si ze< WIOTH  > ) 
height  ASPECT  «  WIDTH) 

get_format I absnaxl ynax  ,  ymin),  nunlen,  fraclen)) 
xl  :*  glxl 0.0  )) 
yg  :■  ymlm 

while  (yg  <■  (ynax  +  0.S  *  yspacing))  do 
begin 

if  abslyg)  <  (0.S  •  yspacing)  then  <•  y  is  rero  •) 

i*  (xmin  <  0,0)  and  ( xmax  >  0.0)  then  <•  x  axis  crosses  y  axis  *> 

yg  :»  yg  +  yspacing)  (•  skip  to  the  next  label  *) 

yl  :*  gly<yg)i  (*  convert  y  from  global  to  local  *> 

movelxl  -  SEMI,  y 1 ) t  (•  move  to  the  proper  location  •> 

line(xJ  +  SEMI,  yl))  («  and  draw  the  tick  nark  *) 
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I’.-'  ■•****  ‘  t 


if  (xmax  <•  0.0)  then  .<•  y  axis  is  at  the  rioht  edge  •> 

move(xl  ♦  WIOTH,  yl  -  0.2S  •  height) 
else  <*  y  axis  somewhere  loft  of  it  •> 

move(xl  -  WIOTH  •  (nunlen  +  2),  yl  -  0.25  «  height  >* 

writenumlyg,  nunlen,  fraclen)i 
yg  !*  yg  ♦  yspacing 
end 

and*  < *  labelyaxis  *  ) 


procedure  plot  (var  x,  y  :  plot„array*  n  :  integer*  linetype  :  integer* 


title  :  plot_string)i 

. . . . * . . . . . . 

<•  *  ) 

<•  Abstract:  This  routine  plots  the  input  arrays  as  a  line  graph.  The  ♦  ) 
(•  x  array  should  contain  the  horizontal  coordinates  and  the  y  array  *) 

<•  should  contain  vertical  coordinates,  n  is  the  total  number  of  • )' 

<•  points,  which  must  be  greater  than  zero,  linetype  is  an  integer  •) 

<«  specifying  the  desired  line  style  -  •) 

< •  t  -  solid  1 lne  • ) 

<  *  2  -  dashed  line  • ) 

( •  7  -  dotted  line  • ) 

<*  < see  page  325  of  the  HP  Pascal  3.0  Graphics  Techniques  Manual  for  •> 

(•  more  styles),  title  specifies  the  title  for  the  plot,  unless  it  •) 

<*  is  declared  a  null  string,  in  which  case  the  current  plot  is  *> 

<•  displayed  over  the  previous  one.  •  ) 

<  •  *  ) 

. . . . . . . . .  ) 


var 

xspaclng,  yspacing  :  real*  <•  spacing  for  tick  narks  •) 

begin 

if  lnitstate  and  (n  >  0)  then 
begin 

if  (title  <>  *')  then  <*  a  new  plot  ») 
begin 

set_l lne_sty le< 1  )  i 
»ri testr< title  )  i 
drawbox ( 1  ,  1  ) » 

aet_v t ewport ( 0 . I  2  ,  0.88,  0,07,  0.65 ) * 
f lndnaxmin! x  ,  n,  xnax  ,  xrun)i 
f indmaxminfy  ,  n,  ynax  ,  ymin)* 
f indspacing( xnax  ,  xmin,  xspacing)* 
f indspacingl ynax  ,  ymin,  yspacing)* 
drauiaxes  i 

labelxaxis( xspacing  )  t 
labelyaxis! yspacing ) 
endi 

set_l ine_atyle( 1 lnetype  )  i 
drawl ines<  x ,  y ,  n  ) 
end 

endi  <  *  plot  *  ) 


end.  <  »  plotpac  *  > 


<*  •> 

<«  USER  OATA  FOR  NEWMARK  INTEGRATION  PROGRAM  *> 

<#  (SIX-DOF  IN-BORE  PROJECTILE  MODEL)  •> 

( •  • ) 

. . . . . . * . . 

. . . . •  ♦  GLOBAL  VARIABLE  DECLARATIONS  . . . . > 

0  '■  real i  <  gravi tationel  acceleration  conetant  > 

mp  i  reali  <  ness  of  projectile  > 

ng  :  reali  •(  mass  of  gun  ) 

Ip  :  reali  <  maos  moment  of  inertia  of  projectile  > 

Ig  i  reali  <  naae  nonant  of  inertia  of  gun  > 

R  :  reali  <  redial  clearance  between  gun  barrel  and  projectile  > 

nu  i  reali  <  coefficient  of  friction  at  bourrelet  (zero)  > 

Lp  i  real i  {  distance  between  obturator  and  bourrelet  > 

Lb  :  real i  <  distance  between  projectile  c.g.  and  bourrelet  > 

Lo  :  reali  1  distance  between  projectile  c,g.  and  obturator  > 

Ka  :  reali  <  angular  stiffness  at  the  obturator  > 

Kb  :  reali  <  bourrelet  stiffness  (linear)  > 

Ko  :  reali  <  nylon  obturator  stiffness  (linear)  > 

Ko!  :  reali  <  aluminum  stiffness  at  obturator  > 

Cl  :  reali  <  angular  damping  coefficient  > 

C2  ■'  reali  <  nonlinear  angular  damping  coefficient  > 

P  :  reali  <  maximum  applied  pressure  to  projectile  > 

area  :  reali  <  cross-sectional  area  of  projectile  ) 

. . .  INPUT  FUNCTIONS  . . . 


function  get_ndof s  »  integer! 

(«  set  the  number  of  dofs"  •> 
oegin 

get_ndofs  <*  8 
ondi  ( *  get_ndof s  *  > 

procedure  get_inj t„values  ( ver  displ  ,  veloc  ,  ac Ctrl  :  vector  )i 

(•  specify  the  initial  values  *) 

var 


yo  J  reel  i 
begin 


zerovl di spl > i 

zerov( veloc )  i 

zerovl accel > I 

g  ■ 

-  386.4i 

< 

in/sec“2  > 

mp 

:•  16.0  /  oi 

< 

lb-sec "2/in 

) 

mg 

I-  4000,0  /  gi 

< 

lb-sec "2/in 

> 

Ip 

0.583i 

{ 

lb-in-*ec'>2 

i 

/ 

lg 

150,01 

{ 

lb-in-sec',2 

/ 

R  i 

-  0.0052i 

{ 

in  ) 

Lp 

! m  5 . 5  i 

< 

in  > 

Ka 

3.6E5i 

< 

in-lb/rad  > 

Ko 

i-  5.0ES) 

{ 

lb/in  > 

C) 

i «  71,0i 

< 

Ib-in-sec  (about  5  percent  of  critical)  ) 

C2 

i  -  0.0i 

P  : 

-  S0.0E3i 

{ 

lb/in'2  > 

area  i-  1 7 . 53 i 

< 

in '2  > 
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writet ’Enter  the  distance  Prom  the  obturator  to  the  c.g.  (in.)  :  '  >i 
readln(Lo  >f 

wrlte( ’Enter  the  initial  projectile  angle  (radians)  :  * >i 
readlnt  diapl£6 3 )» 

write* ’Enter  the  initial  vertical  displacement  at  the  obturator  (in.)  :  *  >  * 
.  reedln( y o >i 
Lb  :*  Lp  -  Loi 

displtS3  »«*  yo  +  Lo  •  sin( disc  1 1 6 3 > i 
end*  (*  get_init_valuea  •) 


procedure  get_maes  (displ, 
<•  specify  the  mass  matrix 
var 

real  i 
real  » 
real  i 
real  i 
:  real i 
real  i 


Xg,  Xgdot 
Yg,  Ygdot 
0,  Odot 
xp ,  xpdot 
yp .  ypdot 
a  .  adot 
begin 

Xg  :*■  displC  1  li 
Yg  displt23 1 

0  !■  displt33i 
xp  displ£43i 

yp  s-  displtS3i 
a  !•  displ£63i 


veloc 

* ) 


J  vector j  ver  mass  >  matrix >i 


Xgdot 

Ygdot 

Odot 

xpdot 

ypdot 

adot 


s*  ve loci  I]  i 
!•  veloc£23i 

•  veloct33i 
J"  veloc£43i 
t"  velocCSIi 

*  veloct63i 


zerom<  mass  )  i 


massC 1,1] 
masst 1  ,3  3 
masst 1,4] 
mass! 1  ,5  3 


mg  +  mp i 

-mp  *  \xp»sin(0>  +  yp«cos<0>)» 
mp  *  cos(0>i 
-mp  •  sin<0)( 


ma  a  s  f  2  , 1 3 
masst 2  ,3  3 
masst  2,43 
masst  2  ,53 


mg  +  mp i 

mp  •  <xp*cos(0)  -  yp«sin<0)/i 
mp  •  sin(0)i 
mp  *  cos(0>i 


masst  3,13 
masst  3  ,2  3 
masst  3  ,3  3 
masst  3  ,4  3 
masst  3  ,5  3 


-mp  *  <xp*sin<0>  +  yp*cos(0>>» 
mp  *  (xp»cos(0)  -  yp»sin(0>>i 
Ig  +  mp  •  (sqr(xp)  ♦  sqr(yp))i 
-mp  *  yp i 
mp  *  xpi 


masst 4 , 1 3 
masst  4,2  3 
masst  4  ,3  3 
masst  4,4] 
masst  4  ,5  3 


mp  i 
mp  i 

mp  *  (xp*(cos(0)  -  sin(Ol)  -  yp*(sin(0)  ♦  cos(0)))i 
mp  *  (sln(0)  ♦  cos(0>)i 
mp  *  ( cojI 0  )  -  s in< 0  > )  i 


masst  5 , 1 3 
masst S  ,2  3 
masst  5  ,3  3 
masst  5  ,5  3 


-mp  *  sir.(0>» 
mp  •  cos(0)i 
mp  *  xpi 
mp  i 


masst  6,6  3  : *  Ip 

endi  (  •  get_mess  •  ) 
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procadura  get_damping  (displ,  veloc  s  vectori  var  damping 
(•  specify  the  damping  matrix  *> 


var 

X g |  Xgdot  J  real i 
-  - Yg ,  Ygdot  :  real  I 
0,  Odot  :  reali 
xp ,  xpdot  :  raali 
yp  ,  ypdot  J  real i 
a,  adot  :  raali 
begin 

rXg  : “  displt 1  1 1 
Yg  :«•  displt23i 
0  :*  displt31i 
xp  :•  displ(4]| 
yp  :•  displtSli 
a  : -  displt 6 ] ( 


Xgdot  valoct 1  ]  i 

Ygdot  :■  veloct21t 
Odot  :■  veloct3)i 
xpdot  veloc[4]| 

ypdot  veloctSli 

adot  :•  velocCSli 


matrix  )i 


zeromt damping )  i 


damping! 1,4]  :•  -2.0  *  «p  *  sin<0>  •  Odoti 
damp l ng t 1,5]  : •  -2 . 0  •  mp  •  cos ( 0 )  •  Odoti 


dampingt  2,3] 


2.0  *  mp  *  < xpdot *coa< 0 >  -  ypdot*ain< 0 > ) l 


dampingt  3,1] 

damping!3 ,2 1 

dampingt 3 ,31 

dampingt  3,4] 
dampingt3 ,51 


i •  -mp  *  ( < xpdot*sin< 0 )  ♦  ypdot*coa( 0 > ) 

+  0dot*( xp*coe< 0  )  -  yp*»in(0>)>i 
mp  *  (( xpdot •cost  0 >  -  ypdot *8in< 0 > ) 

-  0dot*< xp*sin< 0 >  +  yp«cost0>>)i 

-mp  •  ( ( Xgdot »yp  -  Ygdot»xp>  *  ain(O) 

-  <Xgdot*xp  +  Ygdot*yp)  *  costO)>i 
:<•  mp  •  <Xgdot«ain<0>  -  Ygdot«cos(  0  )  > » 

mp  *  ( Ygdot *ain< 0  )  +  Xgdot«coat 0 > > I 


dampingt  4,1] 
dampingt  4,2] 
dampingt  4,3] 
dampingt  4,4] 
dampingt  4  ,5  ] 


mp  »  sin(0)  •  Odoti 

-mp  *  cob(Q)  *  Odoti 

-mp  •  (xp*0dot  +  ypdot)  i 

2.0  •  mp  •  Odot  *  tcos<0)  -  sin(0))i 

-2.0  *  mp  *  Odot  •  (sintO)  *  cos(0>>i 


damping! 5,1) 
damp l ngt  5  ,2  ] 
dampingt  5,3] 
dampingtS ,4] 


-mp  *  coot  0  )  •  Odot i 

-mp  *  s l n( 0  )  •  Odot i 

mp  *  (  Xgdot  "cost  0  >  ♦  Ygdot*sin<  O')  >  i 

2.0  •  mp  *  Odoti 


dampingtG.61  Cl  +  CZ  *  abstadot 
endi  ( *  gat_damping  *) 


procedure  get_st i f fness  (displ,  veloc  ’•  vectori  var  stiffness  :  matrixli 

(*  specify  the  stiffness  matrix  «) 

var 

Xg  ,  Xgdot  '■  real  i 
Yg,  Ygdot  •  reait 
0,  Odot  :  reali 
xp ,  xpdot  :  real i 
yp  ,  ypdot  :  real i 
a  ,  ado t  :  rea 1 i 
Ro  ,  Rb  ••  reel  i 
begin 
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Xg  J»  diapl£13i  Xgdot  »•  velocClli 

Yg  i “  diapl[23i  Ygdot  »■  veloc£23t 

0  :■  displt33i  Odot  «■  veloct33i 
xp  i -  diaplC43i  xpdot  !■  veloc(43i 

yp  J-  dl sp  1 C 5 3 1  ypdot  veloc£53i 

a  »■  di sp 1 1 6 1 f  adot  :•  veloc(E3i 

2erdm( ati f fnaee  > i 

Ro  i"  yp  -  Lo  «  sin(a)« 

Rb  i»  yp  ♦  Lb  •  ain(a  )i 


if  (aba(Ro)  >■  R>  than  (*  contact  at  obturator  «> 

Kol  !-  ) .0E8 

•  lee 

Kot  »■  0.0i 

if  <aba(Rb>  >•  R)  than  (*  contact  at  bourralet  •> 

begin 

Kb  S.0E7i 
nu  :*  0.0 
end 
else 
begin 

Kb  :■  0.0i 
nu  s ■  0. 0 
andi 

ati f fneoaC t  ,43  :■  -np  *  coa<0>  *  aqr(0dot)i 
ati ffneaat 1  ,53  :■  np  •  ain(O)  *  aqr(0dot)i 

ati f fneaa[2 ,43  ! •  -np  •  ain<0>  *  aqr(0dot)i 
ati f fneaaC2 ,53  s*  -np  •  coa(O)  •  eqr(0dot>« 

at i f fneaat 3 ,4 3  :■  2.0  •  np  *  Odot  •  xpdot i 
at i f fneaaC 3 ,53  :■  2.0  *  np  •  Odot  »  ypdoti 

at i f fneaat 4  ,4 1  :■  -np  *  aqr(Odot)  •  <coa(0>  +  ain(0>)i 
at i f fneaat 4  ,5 3  !•  np  •  aqr(Odot)  *  <ain<0)  -  coa(0))i 

at  i  f  fneaat  5 ,5  3  J*  -np  •  aqr(Odot)  +  Kb  +  H'c  +  Koli 
at l f fneaat 5 ,6 3  :■  Kb  *  Lb  -  (Ko  Kol)  *  Lo« 

at i f fneaaC 6  ,5 3  :•  Kb  *  Lb  -  <Ko  +  Kol)  •  Loi 
at i f f neaal 6  ,6 3  :■  Kb  *  aqr(Lb)  +  (Ko  +  Kol  )  *  aqr(Lo)  +  Ka 
endi  (♦  get_at 1 f f neaa  •) 


procedure  get_force  (diapl,  veloc  :  vector*  tine  :  realt'var  force  :  vector >i 

(»  apecify  tha  force  vector  <*  ) 

var 

Xq  ,  Xgdot  :  rea 1 1 
Yg  ,  Ygdot  !  real i 
0,  Odot  :  reali 
xp  ,  xpdot  :  real i 
yp ,  ypdot  :  real i 
a  ,  adot  :  real i 
load,  F  :  real) 

function  agn(x  :  real)  :  reali 
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tarn  the  sign  of  x  •> 


begin 

if  x  '  0.0  then 
egn  s  *  -1.0 

else 

sgn  i m  1.0 
end*  <•  sgn  •) 

begin 

Xg  ; •  displl 1 li 
Yg  !*  displC21i 
0  !•  displt33« 
xp  J-  displtili 
yp  !■  displt 51 i 
a  !•  displCSh 

zerov( force  >i 


Xgdot  *■  velocC 1  ] i 
Ygdot  •*  velocC23i 
Odot  }■  velocl3]» 
xpdot  «•  v*iooC41i 
ypdot  »•  veiocC5]i 
adot  »”  veloe[B)i 


if  (tine  <  0.0025)  then 

load  -  0.1  +  500.0  *  tine 
else 

load  >•  1.5  -  180.0  •  tine t 


i  f  (  load  >  1.0)  then 
load  i*  1.0 
else 

if  (load  <  0.0)  then 
load  s*  0.0i 

F  :•  load  •  P  •  area i 

force! 21  :•  -g  •  <np  +  ng)| 

forcaC4]  -np  *  g  *  sin(0)  -  nu  •  Kb  *  (abs(yp  *  Lb»a )  -  R)  *  sgn(xpdot) 
♦  F  •  coe(a )» 

Porcel5]  i-  Kb  •  R  •  sgn( yp  ♦  Lb  •  a )  -  np  •  g  »  coa(0>  +  F  •  sln(a> 

+  R  *  Kol  «  sgn(yp  -  Lo  •  a>» 

forcetS)  :•  Kb  *  R  *  Lb  •  ogn(yp  +  Lb  •  a) 

-  R  •  Kol  •  Lo  •  egn<yp  -  Lo  *  a  > ( 

endi  ( •  get_force  *  > 
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<•  USER  DATA  FOR  NEWMARK  INTEGRATION  PROGRAM  •> 
(  »  *  ) 
. . . . .  ) 


x ••••*••••••••••••••••••••••••  INPUT  functions  . . . . . 

function  get_ndofs  s  integer) 

<•  set  the  number  of  dofs  •) 
begin 

get_ndof»  :■  2 
-""endi  (•  get_ndofs  •) 


procedure  get_init_values  (var  dispi,  veloc  ,  accel  :  vector)! 

<•  specify  the  initial  values  *) 

begin 

zerov( displ>i 
zerov( veloc ) i 
zerovt  accel ) i 
accel£2)  :■  10.0 
endi  (*  get_init_values  •) 


procedure  get_mass  (dispi,  veloc  :  vectori  var  mass  :  matrix)! 

<•  specify  the  mass  matrix  •) 

begin 

zerom<  mass )  i 
mass! 1,13  :•  2.0i 
mess[2,23  :•  1.0 
end!  (•  gat_mass  «> 


procedure  get_damping  (dispi,  veloc  s  vector!  var  damping  :  matrix)! 

<»  specify  the  damping  matrix  *) 

begin 

zerom<  damping ) 
end!  <*  got_damping  •) 


procedure  get_stif fness  (dispi,  veloc  :  vector)  var  stiffness  :  matrix)) 

<*  specify  the  stiffness  matrix  *) 

begin 

zerom( stiffness ) t 

stif fnesstl  ,1 1  :■  6.0)  st l f fnessC 1  ,2 3  :«-2.0i 
st if fnesst2  , )  3  :*-2.0i  st l f fness(2  ,2 3  :■  4.0 
end)  (•  get_st l f fness  *) 


procedure  get_force  (dispi,  veloc  :  vector)  time  '■  real)  var  force  :  vector)) 

(*  specify  the  force  vector  *) 

begin 

zerov( force  )  i 
forcaI23  :■  )0.0 
end)  (*  get_force  *) 
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LIST  OF  SYMBOLS 


an  coefficient  of  nth  term  In  power  series  representation  of  the 

nonlinear  foundation  moment 

Ap  area  of  projectile  cross  section  on  which  gas  pressure  acts 

Cj,  c2  first  and  second  order  damping  coefficients  of  projectile  resistance 
to  yawing  motion 

fC]  damping  matrix  In  system  equations  of  motion 

f  ,  f  applied  forces  on  projectile  In  x  and  y  directions 

x  y 

Ff  friction  force  on  projectile  at  contact  between  bourrelet  and  gun- 

Tbc  bore 

(F)  column  vector  of  forces  In  system  equations  of  motion 

g  acceleration  due  to  gravity;  also  denotes  gun 

G  denotes  gun  CG 

1,  j  unit  vectors  In  X,  Y  directions 

Ig  mass  moment  of  Inertia  of  gun  about  an  axis  through  Its  CG  perpen¬ 

dicular  to  plane  of  motion 

Ip  mass  moment  of  Inertia  of  projectile  about  an  axis  through  Its  CG 

perpendicular  to  plane  to  motion 

kpc  spring  constant  for  bourre let/gun-bore  contact  spring  deflection 

model 

kQ  stiffness  of  obturator  plastic  band  In  obturator  transverse  spring 

model 

I 

kQ  stiffness  of  metallic  part  In  obturator  transverse  spring  model 

[K]  stiffness  matrix  In  system  equations  of  motion 

Instantaneous  distance  of  bourrelet  from  gun  CG  measured  along  gun 
bore  axis 

tp  distance  between  projectile  CG  and  plane  of  bourrelet 

iQ  distance  between  projectile  CG  and  plane  of  obturator 

i,  m  unit  vectors  along  x#  y  directions 


77 


L 


m(t) 

CM] 


P 

p(t) 

P 

6P 

(q) 


Rcl,  be 

Rcl,o 

t 

T 


Vk 

V 

x,  y,  z 


X,  Y,  Z 


Lagranglan  function 

mass  of  gun 

mass  of  projectile 

applied  yawing  moment  on  projectile 

Inertia  matrix  In  system  equations  of  motion 

number  of  terms  In  power  series  used  to  represent  the  nonlinear 
foundation  moment 

denotes  projectile 

time  history  of  gas  pressure  on  projectile 

denotes  projectile  CG 
virtual  power 

rth  generalized  coordinate 

column  vector  of  generalized  coordinates  In  system  equations  of 
motion 

nonconservative  generalized  force  associated  with  rth  generalized 
coordinate  qr 

radial  displacement  vector  of  projectile  centerline  In  plane  of 
bourrelet  measured  from  the  centerline  of  the  gun  bore 

radial  clearance  between  bourrelet  and  gun  bore 

radial  clearance  between  obturator  and  gun  bore 

time 

kinetic  energy 

generalized  velocity  In  direction  of  kth  generalized  coordinate 
potential  energy 

gun  body-fixed  axes  of  coordinates 

coordinates  of  projectile  CG  relative  to  gun  body-fixed  axes  of 
coordinates 

Inertial  frame  of  reference 
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Xg,  Yg  displacement  of  gun  C6  relative  to  Inertial  frame  of  reference  along 
X,  Y  axes 

a  angular  (yawing)  displacement  of  projectile 

«ljC  displacement  of  projectile  into  the  gun  bore  at  contact  between 
bourrelet  and  gun  bore 

60  projectile  displacement  at  obturator 

e  angular  displacement  of  gun  In  direction  of  positive  Z  axis 

yb  coefficient  of  friction  for  bourrelet/gun-bore  contact 
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